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Abstract 

We present a general approximation framework for computation of optimal feed- 
back controls in linear quadratic regulator problems for nonautonomous parabolic dis- 
tributed parameter systems. This is done in the context of a theoretical framework 
using general evolution systems in infinite dimensional Hilbert spaces. We discuss 
conditions for preservation under approximation of stabilizability and detectability 
hypotheses on the infinite dimensional system. The special case of periodic systems is 
also treated. 
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1 Introduction 


In this paper we present a theoretical approximation framework for computation of op- 
timal feedback controls in linear quadratic regulator problems governed by parabolic 
partial differential equations with time dependent coefficients. Our efforts were orig- 
inally motivated by the desire to develop control strategies (distributed in nature) 
for insect dispersal models (see Chapter 1 of [BK2] and the references therein) which 
have been shown to involve time dependent coefficients. 

The presentation below is somewhat in the spirit of that for autonomous 
parabolic systems in [BK1] and [LT] in that we attempt to develop a convergence 
theory in which uniform stabilizability of the original system is preserved under ap- 
proximation. It differs substantially from [BKl] and [LT] since we don’t use directly 
sectorial properties of the operators and resolvent and spectral set arguments to es- 
tablish preservation of stabilizability and detectability. (Indeed, the time dependent 
nature of our system prevents this.) Nor do we use the Trotter-Kato theorem (which 
is not well suited for use with nonautonomous control systems) in our convergence 
arguments. 

In section 2 we summarize previous results for abstract LQR problems on 
infinite time intervals and formulate these in a form readily used in our subsequent 
discussions. This formulation is based on the abstract frameworks found in [CP], [G], 
[BKl] and [Da], [Dll], [DI2]; we rely heavily on the ideas of Da Prato and Ichikawa 
which guarantee uniqueness of solutions of the associated Riccati integral equations 
under certain stabilizability and detectability assumptions. 

An approximation framework for abstract evolution systems in the spirit of 
[G], [BKl] is given in section 3; convergence of the approximate Riccati operators 
( and, of course, the corresponding controls and trajectories) is established under 
uniform stabilizability and detectability hypotheses on the approximate evolution 
control systems. 
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Our major contributions are given in section 4, along with the presentations 
in section 5 and 6, where we show how the hypotheses of section 4 can be verified 
for rather wide classes of problems of interest. In section 4 we focus our attention 
on parabolic systems described by time dependent sesquilinear forms (in the spirit 
of the autonomous system frameworks in [BK1], [BI1], [BI2]) and associated evolu- 
tion equations. We make substantial use of the results of Tanabe [T] to formulate 
our problems in a weak (V*) sense. Our fundamental convergence results (Theo- 
rem 4.4) for the uncontrolled systems rely on a sesquilinear or variational formulation 
of the systems, strong V -ellipticity of the parabolic evolution systems, approximation 
properties for the spaces approximating the state space, and the Gronwall inequality. 
(Certain aspects of this approach can be relaxed to allow us to treat weakly damped 
hyperbolic systems-see [BKS], [BKW].) We are then able to reduce convergence 
questions for the controlled systems (e.g., convergence of Riccati variables, optimal 
controls and feedback evolution systems) to conditions of uniform stabilizability and 
uniform detectability of the approximate systems (Theorem 4.5). 

We show in section 5 that we can obtain these uniform stabilizability /detectability 
conditions by preservation under approximation of dissipative inequalities for certain 
classes of evolution control systems. Sufficient conditions that are readily checked in 
many examples are given and several special cases are noted. 

An alternative approach is presented in section 6 where we restrict our consid- 
erations to parabolic systems for which the domain V of the generator of the evolution 
system embeds compactly in the state space H. In this case, it is shown that stabi- 
lizability / detectability of the original system is preserved under approximation. 

Finally in section 7 we give an example of a class of parabolic partial differential 
equation control problems for which all the hypotheses of our theoretical framework 
can be easily verified. 

We have used the ideas presented in this paper to develop and test compu- 
tational packages for solving nonautonomous parabolic control problems of the type 
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discussed in section 7. However, since our presentation here is already quite long and 
since a presentation of our detailed numerical findings would entail lengthy discus- 
sions, we will not discuss the numerical examples. A separate manuscript is under 
preparation; the interested reader can also consult [W]. 

We believe that the present paper offers new results for time dependent infinite 
dimensional control systems. Moreover, our arguments are such that we offer an 
attractive alternative approach to those found in [G], [BK1], [LT] even in the case of 
autonomous parabolic systems. 

2 The abstract linear quadratic regulator problem on an 
infinite time interval 

In this section we formulate a linear quadratic regulator problem for evolution system 
dynamics in a Hilbert space. We present a collection of functional analytic and 
control theoretic results related to such problems. The results we give in this section 
are known and, while in some cases we have modified the statements to present the 
results in a form most suited to our purposes, the reader can easily refer to the 
literature for proofs. In particular, we use freely results found in [CP] and [G] and 
rely heavily on recent results of Da Prato and Ichikawa [Da], [Dll], [DI2]. 

We first recall results for evolutionary systems. Let H be a Hilbert space 
with inner product <•,•>. Let A(t 0 ,tj) = {(L-s) | to < s < t < t/}, A oo (i 0 ) = 
{(t,s) \ t 0 < s < t < 00} and L(H ) be the Banach algebra of bounded linear operators 
on H. We use -^([fo, */]; L(H)) to denote the set of operator valued functions that 
are bounded on [to,t/]. We recall that T(-,*) : A (to,tf) L(H) is called an evolution 
operator if T satisfies the following conditions: (i) T(t, s) = T(t , r)T(r, s), for to < s < 
r < t < tf] (ii) T(t,<) = /, for t 6 [to,tj]-, and (iii) T(t,s) is strongly continuous in s 
on [< 0 , t] and strongly continuous in t on [s,f/]. We say that an evolution operator has 
exponential growth if there exists Mi > l,u > 0 such that ||T , (f,s)x|| < Afie^ 1- ^||x||, 
for (t,s) € A 00 (^ 0 ), # € H. An evolution operator is said to be uniformly exponentially 
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stable if there exists M > 1 and a > 0 such that ||T(f, .s):r|| < Me -0 '( <-4 )||x||, for 
( t,s ) G Aoo(< 0 ),a: G H. We have the following fundamental results of Datko which 
will be crucial to our presentation. 

Lemma 2.1 [Dt] Consider an evolution operator T { •,•) with exponential growth. 
Then is uniformly exponentially stable if and only if there exists an M 2 such 

that 

/ OO 

\\T(t, s)x\\ 2 dt < M 2 ||x|| 2 , for s >t 0 ,x£ H. 

Furthermore we can find constants M3 > l,a > 0 depending only on A/ 1? M 2 and u? 
for which the following estimate holds: 

\\T(t, 5)|| < M 3 e~ a ^~ s \ for (t, s) € AM. 


The original statement and proof of this theorem are due to R.Datko. We 
have modified slightly (see Appendix A of [W]) the original proof in [Dt] to point 
out the relationship between the constants A/3, a and A/j,A/ 2 and u>. This will be 
essential for our subsequent use with approximation systems. In [W] it is shown that 
the constants M3 and a can be chosen as: 


A/ 3 = 2 A/j e 4M * Mi M2wAf a +i) ^ 


log2 

AM 2 M 1 2 ( 2 ojM 2 + 1 )' 


In our discussions of control systems, perturbations of evolution operators 
(see [CP]) will play an important role. Let tj < 00. Consider a uniformly bounded 
evolution operator and C(-) G 5oo([*0i / /] ; L(H)). Then the integral equation 

for 5(f,5) G L(H ) given by 

S(t,s)x = T(t,s)x + J T(t,T))C(r))S(r),s)xdr], for x G H, 

has a unique solution 5 (-,-) in the class of strongly continuous operator valued func- 
tions. Moreover, S( •,•) is an evolution operator and is called the perturbed evolution 
operator corresponding to the perturbation of T(-,-) by C(-). In addition, 5 '(-,-) is 
also the unique solution of 

S(t,s)x = T(t,s)x + / S(t,Tj)C(T])T(Tj, s)xdrj, for x G H. 
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We turn next to our formulation of the regulator problem for an evolution 
system. We let H,U be real Hilbert spaces with inner products < >//, < >u\ 

H , U are the state space and the control space respectively. Consider an evolution 
operator T(-,-) defined on A oo (^ 0 )- The control system is described as follows: for 
any u 6 L 2 ([t 0 , oo); U ), the corresponding trajectories satisfy the equation: 

(2.1) x(t) = T(t,s)x(s) + J T(t,T)B(r)u(r)dT , for ( t,s ) 6 (t 0 ). 

The cost functional is given by 

r OO 

(2.2) Joo(u; to, xo) = / {< W(t)x(t),x(f) ># + < R(t)u(t),u(t) > v } dt, 

Jt 0 

where a:(-) is the trajectory corresponding to u with x(t 0 ) = x 0 . For each given t 0 ,x 0 , 
the optimal control problem is to find a control u* which minimizes (2.2) over all 
u € L 2 ([to, oo); U). 

We can consider (2.2) as the limit as tf. — > oo of 

(2.3) J(u; t 0 , t k , x 0 ) =< Gx(t k ), x(t k ) > H 

+ / {< W(t)x(t),x(t) >H + < R{t)u(t),u(t) >u}dt, 

Jto 

with G — 0. Here we shall summarize existence results for optimal controls in the 
infinite time interval, existence and uniqueness of the solutions of the Riccati integral 
equation on an infinite time interval, and stability of the feedback system. 

We make the following standing assumptions for all subsequent discussions 
of (2.1), (2.2): (i) The evolution operator T(-,-) has exponential growth. (Thus, in 
particular, T(t, s ) is uniformly bounded for s, t in any bounded sub-interval of [f 0 , oo)); 
(ii) The operator valued function B(-) : [fo>oo) i-+ L(U,H ) is uniformly bounded in 
[f 0 , oo), i.e. there exists Mb such that \\B(t)\\mj ik f) < Mg for all t € [£o> oo); (iii) The 
operator valued function W(-) : [< 0 , oo) >-*• L(H ) is uniformly bounded in the interval 
[< 0 , oo), and W(t) is nonnegative definite self-adjoint for all t 6 [to ? °o); (iv) The 
operator valued function R(-) : [< 0 , oo) i— ► L(U) is uniformly bounded in the interval 
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[t 0 , oo), and R(t ) is positive definite self-adjoint for all t £ [t 0 , oo). Furthermore, there 
exists a constant r > 0 such that < R(t)u,u >(/> r||u||y, for all u £ U. 

Under these assumptions we consider the linear quadratic control problem in 
the interval [t 0 , tk] for tk < oo. That is, we consider the cost functional (2.3) with our 
system (2.1). Then for any bounded self-adjoint nonnegative definite linear operator 
G, it is well known that for each given xo £ H , there exists a unique control u such 
that 

J(u-,t 0 ,tk,x 0 ) < min J(v,t 0y tk,xo). 

This control u can be written in a feedback form u(t) = — R' 1 (t)B*(t)Q(t)x(t), for 
t G where i(-) is the corresponding trajectory and Q(-) : [to,tf\ t-> L(H), is 

the unique self-adjoint solution of the Riccati integral equation (RIE) 

(2.4) Q(t)x = T m (t k ,t)GT(t k ,t)x + jf T*(r],t)W(ri)T(r),t)xdT) 

~ J t T*(ri, t)Q(ri)B(ri)R~ 1 (T])B*(T])Q(T])T('r], t)xdr) 
for all t G [to, tfc] and x G H. 

We note that in the case G = 0, the above equation reduces to 

(2.5) Q(t)x = j‘“ r(q,t)[WW - Q(n)Bh)R-'(<,)B-WQ(v))T( v ,t)xd n , 

for all t G [to, h] and x G H. We also note that the equation (2.4) is equivalent to 

(2.6) Q(s)x = T*(t,s)Q(t)T(t,s)x + J T*{r),s)W(r))T(ri,s)xdr] 

- J‘ T ‘ (,, s)Q( n )B( v )R-'(i)B-( V )Q(n)T( v , s)xi v 
Q(tk)x = Gx 

for all t 0 < s < t < tk and x G H. Solutions of this latter equation have a rep- 
resentation that is often used in control theoretic arguments. Consider any u(-) € 
T 2 ([t 0 , tfc]; £/), and for x G H, define a .//-valued function y(- ) by 

y{t) = T(t,s)x + J T(t,T)B(T)u(T)dr, for t G [5,t fc ]. 


6 


If Q(- ) is a self-adjoint solution of (2.6), then 
<Q(s)x,x> h - < Gy(t k ),y(t k ) > H 

(2.7) + jf ‘ {< W(t)y(t), y(t ) > H + < R{t)u(t), u(t) >t/} dt 

— j < R(t)z(t),z(t) >u dt, 

where z(t) = u(t) + R~ 1 (t)B*(t)Q(t)y(t). This can be used to show that (2.6) has a 
unique self-adjoint solution. 

Before continuing our discussion, let us introduce additional notation. Let 

£+ = | E | E 6 L(H),E self-adjoint, nonnegative definite. } 

C s ([t 0 , <*]; E + ) = | K : [<o>4] •— 1 * S + | K strongly continuous, j 


By the uniqueness of the solution of equation (2.6), we can define a mapping A : 
£ + C s ([<o, <*]; S + ) as following: for each G € £ + , A G is the unique nonnegative 

definite self-adjoint solution of equation (2.6). Under our general assumptions, it is 
easily seen that for fixed G the map A depends only on f*; if we consider the linear 
quadratic regulator problem on two bounded intervals [to 5 U] and we will use 

Aj, A 2 to denote the maps associated with each interval respectively. 

Now consider a increasing sequence 1? with tk < 00. The map A*, as- 

sociates with each finite interval problem the Riccati equation on [to,tk]. Let G = 0 
and Qk(-) = AfcCr. For simplicity, consider a bounded interval [a, b] C and for 

each t 6 [a, b],x G H , we assume that there exists a constant M(t,x) such that for 
all k 

(2.8) < Qk(t)x,x >h < M(t,x). 

The following theorem (see [Da], [G]) establishes the connection between control 
problems on a finite time interval and problems on an infinite time interval. 

Theorem 2.1 Under the above assumptions, we can conclude: 
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(i) For each t E [a, 6], there exists a unique operator Q(t ) E E + such that Qk(t) — ► 
Q(t) strongly and the convergence is uniform in [a, 6]. Therefore Q(-) is strongly 
continuous, so uniformly bounded in [ a,b ]. 

(ii) y 4 s a consequence of (i), we can define the perturbed evolution systems Sk(-, •)> $('■> •) 
corresponding to the perturbation ofT ( •,•) by —BR~ 1 B*Qk and —BR~ l B*Q re- 
spectively. We have Sk{t,s)x — > S(t,s)x, for all x E H, and a < s < t < b. 
Furthermore the convergence is uniform in t for t E [5,6]. If T( •,•) is jointly 

strongly continuous , then the convergence is uniform for all a < s < t < b. 

1 

The only assumption on the sequence {<*} is that tk increase as a function of 
k. In particular, the above theorem is valid when tk — * 00, as k — > 00. Paralleling 
the usual approach to finite dimensional regulator problems, we can use these results 
to establish results for the control problem on an infinite time interval. To that end, 
consider a sequence with t k — ► 00 as k — * 00. Let (?*(•), 5 jt(-, •) be defined 

as above. If for each t > t 0 we can find a constant M(t) such that < Qk(t)x,x >//< 
M(t)||a;|| 2 , then by Theorem 2 . 1 , we have ( 5 (-) 5 *S'(-, •) defined on [t 0 , 00). Furthermore 
for any (t,s) E A oo (t 0 )» Q satisfies 

( 2 . 9 ) Q(5)a; = T*(t,s)Q(t)T(t,s)x + T*(q,s)W(q)T(ri,s)xdr) 

- f r( v , s )Q( v )BWR-'WB-(v)Q(n)T( v ,s)xd v . 

The equation ( 2 . 9 ) is called the Riccati integral equation (RIE) for the infinite time 
interval. We know from Theorem 2.1 that Q is strongly continuous and uniformly 
bounded in any bounded interval, but Q is not necessarily uniformly bounded in the 
\ entire interval [< 0 , 00). If Q is not uniformly bounded, that implies the minimal cost 

for some initial state x will tend to infinity as tk tends to infinity; that is, there is 
no control yielding finite cost for the infinite time interval problem. Let us state a 
condition which prohibits this situation. 

Definition 2.1 (W-stabilizability) We say that ( 2 . 1 ), ( 2 . 2 ) is W-stabilizable if 
> there exists a constant M such that for any s > t 0 and x E H , we can find a control 


8 


it G L 2 ([t 0 , oo); U) satisfying 

(2.10) Joo(u;s,x) < M||x|| 2 . 

One can then prove (see [Dll], Theorem 3.1) that Q = is a uniformly bounded 

solution of the Riccati integral equation (2.9) in [to, oo) if and only if (2.1), (2.2) is 
W-stabilizable. In this case we have Q(t ) < M • I for t € [to, oo). Furthermore, 
if Q is any other bounded self-adjoint solution of (2.9), we have that Q(t) < Q(t ) 
for t G [t 0 , oo). It follows that using any sequence {tk} with tk —> oo, in the above 
limiting procedure yields the same solution Q to (2.9), which we shall refer to as the 
minimal bounded nonnegative self-adjoint solution of the Riccati integral equation 
on [t 0 , oo) and denote by Q m ,„. 

We note that if the system (2.1), (2.2) is W-stabilizable, then for any s > t 0 
and x G H, the unique optimal control for the infinite time interval problem is given 
by u(t ) = —R~ 1 (t)B*(t)Q(t)S(t,s)x. 

Next we consider a uniformly bounded solution Q of (2.9) and let S be the 
evolution operator corresponding to the perturbation of T by —BR~ l B*Q. We say 

A A 

that Q is a stability solution of (2.9) if S(t,s)x — > 0 as t — > oo for all s > to, x 6 H. 

It is shown in [Dll] that there is at most one stability solution of (2.9). More- 
over, if Q is a stability solution satisfying Q(t) < M • I and Qk is the solution on 
[to, 4] with Qk(tk ) = M • I, then Q(t) < Qk(t) for t G [to, tk] an d Qk(t)x — * Q(t)x as 
k —* oo for each x G H. In addition, if Q is any uniformly bounded solution, then 
Q(t) < Q(t), t G [t 0 , oo); that is, any stability solution is the maximal (uniformly 
bounded) solution. Finally, if the system (2.1), (2.2) is W-stabilizable and if the 
minimal solution Q m i n of (2.9) is a stability solution, then it is the unique uniformly 
bounded solution of the RIE (2.9). 

From the above remarks, it is clear that it is desirable to have conditions that 
guarantee a solution of the RIE be a stability solution. One such condition is a 
detectability condition. 
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Definition 2.2 (W-detectability) Let V(t) = yjw\t). We say that the system 
(2.1), (2.2) is W-detectable if there exists a uniformly bounded function K(-) with 
K(t) € L(H ) such that the evolution operator Tkv corresponding to the perturbation 
of T by KV is uniformly exponentially stable. 

We then have the following result. 

Theorem 2.2 Suppose that the system (2.1), (2.2) is W-stabilizable and W-detectable. 
Then the minimal solution Q m in °f the RIE is the unique uniformly bounded solution 
of (2.9) and the evolution operator S defined by perturbation ofT by — B R~ l B* Q min 
is uniformly exponentially stable. In fact, 

\\S(t,s)\\<Me-^- a \ (t,s)eA oo (t 0 ), 

where the constants M and a depend only on the bounds for B, K, R~* , Qmin and 
MkViP in the bound \\Tfcv(t, 5)11 5: Mkv exp {— fd(t — s)}. 

The first part of this theorem follows from [Dll] (Prop. 3.3). That the con- 
stants M and a depend only on the bounds indicated follows from use of the modified 
Datko lemma, Lemma 2.1 above. As we shall see in the next section on approxima- 
tion, this dependence (or lack thereof) will allow us to infer a uniform exponential 
stability of the approximate feedback control systems whenever we have a uniform 
W-detectability condition satisfied by the approximate systems. 

To conclude this section, we recall that an evolution operator is said to be 
^-periodic if for any (t,s) E A co (< 0 ), we have T(t + 6 , s + 6)x — T(t,s)x , for all 
x E H. We note that any ^-periodic evolution operator satisfies the exponential 
growth assumption that is part of our standing assumptions in this paper. It is 
also easily argued that if the linear quadratic regulator problem is 0-periodic (i.e., 
B,W,R and T of (2.1), (2.2) are 0-periodic), then the minimal solution and the 
stability solution of the RIE are 0-periodic. Of course, we cannot argue that every 
uniformly bounded solution of the RIE is periodic under a periodicity assumption on 
the problem. 
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We turn next to approximation results for the abstract linear regulator problem 
on an infinite time interval. 

3 Approximation of linear quadratic regulator problems on 
an infinite time interval 

Let H n and U N be families of finite dimensional subspaces of the original state space 
and control space H, U respectively. For each N an approximate control system is 
described by 

(3.1) x N (t) = T N (t,s)x N (s) + for (t,s) <E A^o), 

where T N ( •,•) : Aoo(to) t— > L(H N ) is an evolution operator, and B N (-) : U N ► H N . 
The cost functional is given by 

roo 

(3.2) j£(u N ;t 0 ,x%) = / < W N (t)x N (t),x N (t) > H dt 

Jto 

roo 

+ < R N (t)u N (t),u N (t) >u dt 

Jto 

where x N (-) satisfies (3.1) and x N (t 0 ) = Xq . Suppose that each of the approximate 
systems satisfies the standing assumptions for (2.1), (2.2) given above and that each 
is W-stabilizable. Then we can guarantee existence of Q N ( •), the minimal uniformly 
bounded solution of the associated Riccati integral equation on the infinite time in- 
terval [f 0 , oo). Let S N (-, •) be the perturbed evolution operator corresponding to the 
perturbation of T N by — B N (R N ) -1 B* N Q N . In this section, we present results on the 
convergence of Q N , S N . 

We need to make some basic assumptions on the approximate systems. Let 
{H N }x =1 , {U n }n =1 be subspaces of H,U respectively, and P£j,Py be projection 
operators which are assumed to satisfy \\Pfj x ~ x lltf -*■ II Pu u — «||t/ — * Oi as 
N —* oo, for all x G H, u G U. 

We note that the usual orthogonal projections of H and U onto H N , U N re- 
spectively satisfy these assumptions if H N , U N approximate H and U in an appropri- 
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ate sense. (We shall specify approximation systems that satisfy these conditions in 
subsequent sections.) We make the further assumptions on our approximate systems. 

Hypothesis 3.1 (Uniform boundedness) 

(i) There exist constants M > 1 and u > 0 such that 

lir((,s)|U,H) < ||r w (i,s)|| i(H ») < «<■ (■-> 

hold for all N and ( t,s ) E A^to); 

(ii) There exists a constant Kb such that 

\\B(t)\\L(u t H) < Kb, ||5 ;v (<)|| L ( t /N )W N) < Kb 

for all N and t 6 [t 0 , oo); 

(iii) There exists a constant Kw such that 

l|W(t)||i,(//) < K\ v, ||W /Ar (2)||i,(tf*) < Kw 

for all N and t E [to, oo). Furthermore W(t),W Ar (t) are nonnegative definite 
self-adjoint for all t E [to, oo). 

(iv) There exists a constant Kb such that 

||-R(0||l([/) ^ Kri ||-^ Ar (0IU(u w ) — Kr 

for all N and t E [to, oo). In addition, R(t), R N (t) are positive definite self- 
adjoint for all t E [to, oo). There exists a constant r > 0 such that R(t) > r • I, 
R N (t) > r • I, for all t E [to, oo). 

Hypothesis 3.2 (Pointwise convergence) The operators T N (t,s), T’‘ N (t,s), B N (t), 
B* N (t), G N , W N (t), R N (t ) converge strongly to T(t,s), T*(t,s ) , B(t), B*(t), G, 
W(t), R(t ) for any t 0 < s < t < oo, where G,G N are nonnegative self-adjoint opera- 
tors in L(H), L(H N ) respectively. 
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From arguments in [Dll] and [W], it is readily seen that W-stabilizability (i.e. 
condition (2.10)) is equivalent to the following: there exists a constant M > 0, and a 
uniformly bounded feedback operator K(-) : [t 0 , °o) L(H,U ) such that if Tk(-,-) 

is the perturbed evolution operator corresponding to the perturbation of T by BK , 
then for any s > t 0 , x € f/, the cost of the feedback control u(t ) = B(t)K(t)Tx(t,s)x 
satisfies J^u^Syx) < M||a:|| 2 . 

To guarantee the existence of uniformly bounded solutions of the Riccati in- 
tegral equation on the infinite time interval for each of the approximate systems, we 
make a uniform W-stabilizability assumption. 

Hypothesis 3.3 (Uniform W-stabilizability) There exists a constant M > 0 
such that for all N , there exist uniformly bounded feedback operators K N (-) : [^o? oo) 
L(H N ,U N ) satisfying the following: for all s > to and x N (= H N , the feedback control 
u N (t) = B N (t)K N (t)T^- N (t,s)x N has a cost satisfying 

J»(u n -,s,x n )<M mi 2 ,. 

Now consider with — i ► oo as k — ► oo. For each N, let E +Ar be 

the set of nonnegative self-adjoint linear operators in H N ; we define the map : 
E +yv h- f C s ([f 0 , <*]; E +;v ) via the finite dimensional Riccati integral equation on [^ 0 , tk] 
as before. Let G 6 E+, and G N € E +N . Define Qjf(-) = A jfG N , and Q k (-) = 
A k G. Let the evolution operators Sjf and S k correspond to the perturbation of 
T N , T by —B N (R N )~ 1 B* N Q k and —BR~ 1 B*Q k respectively. The theories of the 
approximation of linear quadratic control problems on a finite time interval (e.g. 
see [G], [BK1]) guarantee that under Hypothesis 3.1, 3.2, for each fc, Q k (t ) and 
S k (t,s) converge strongly to Q k (t ) and S k (t,s) respectively as A -> oo for every 
< s < t < t k . Furthermore the convergence is uniform in the interval [<o, t k ], if we 
replace Hypothesis 3.2 by the following assumptions. 

Hypothesis 3.4 (Continuity and uniform convergence) The operator valued 
functions B(t ), B*(t), W(t), R(t) are strongly piecewise continuous in t ( with only a 
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finite number of discontinuity points in any bounded interval); the evolution operators 
T,T* are jointly strongly continuous. The convergences in Hypothesis 3.2 are uniform 
in t and ( t,s ) on any bounded interval. 

From the theory of the linear quadratic control problem for the infinite time 
intervals, there are two cases where Q % converges strongly as k — > oo. In the first 
case, let Qm.n(0 b e the minimal uniformly bounded solution of the Riccati integral 
equation in H N on the infinite time interval [t 0 , oo). If G = G N = 0, then for any N 
we have Q(f ( t)x N — > Q^ in (t)x N , as k — ► oo. Furthermore the convergence is uniform 
in t for t in any bounded interval [t 0 ,tj]. In the second case, we assume the following 
conditions hold: 

Hypothesis 3.5 Assume that there exists a stability solution Q s (-) of the Riccati 
integral equation for the infinite time interval infinite dimensional system and there 
exists a stability solution Q^(-) of the Riccati integral equation on the infinite time 
interval for each approximate system. (Then the evolution operator S N corresponding 
to the perturbation ofT N by —B N (R N )~ i B mN Q ! f satisfies S N (t,s)x — » 0, as t — * oo, 
for s > t 0 , x £ H N .) Furthermore assume that there exists a constant M such that 
for each N , Q^ 1 (t) < M • I for all t > t 0 . Also assume Q s (t) < M • I for all t > t 0 . 

Assuming that Hypothesis 3.5 holds, we let G = G N — M ■ I and , Qk be 
the solutions of the RIE on [t 0 , tk] satisfying Qk(h) = G N , Qk{tk) = G. Then from 
our results for stability solutions given in section 2 we have Qjf ( t)x N — ► (t)x N , 

Q k (t)x — ► Q s (t)x as k — ► oo for all x N € H N and iGF. 

We note that if S N (t, s)x — > 0 uniformly in N, then we have and Q(f (s) 

uniformly bounded for all k and N. To see this, we consider Qk , Q, as given above. 
Then we have, using the relationship in (2.7) with y(t) — S(t,s)x where S is the 
evolution operator corresponding to the perturbation of T by —BR~ l B*Q S , 

< (Qk{s) — Qs{s)) x ,x >h < < {Qk(h) — Qs{tk))y{ik),y{tk) >h 

< 2M\\y(t k )\\ 2 . 
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Since y(tk ) — ► 0, it follows using the uniform boundedness principle that Qk(s) is 
uniformly bounded for all k. Repeating this argument with Q % ( s ), Q ^ (s) and y N (t ) = 
S N (t,s)x, we see that the uniform (in N ) decay of S N yields the claimed uniform 
boundedness for (s). 

In each of the two cases above, we have the following situation: 

QSMPg* Q7(*)P8* 

N -» oo | [7 N oo 

Qk{t)x k ^° Q s (t)x. 

It is desirable in computations to work directly with Q ^ and hence we seek 
results which will guarantee the convergence Q ^ — * Q s of this diagram. To obtain 
such a result, we shall make use of a uniform decay rate for the S N defined via Q 1 ^ . 

Theorem 3.1 Assume that Hypothesis 3. 1-3.3, 3.5 hold. Further assume that for all 
s > to, x € H and e> 0, we can find t such that for all t>t, we have ||5(t, s)x|| < e 
and ||5' ;v (t,s)P^x|| < e for all N. Then Q^^Pffx — ♦ Q a (t)x for all to < t < oo. 

Proof: Let M be the bound for Q $ and Q ? that are the stability solutions of 
Hypothesis 3.5. Let and Qk be the related RIE solutions on [^o? h] satisfying 
( tk ) = M • 7, Qkifk) = M • I. Then for t <tk we have 

II Q?(t)PKx - Q.(t)x II 2 < ||(<?f(0 - Ql(t))PK* II 2 + WQHWPhx - 

+IIW.W - <W())*II 2 - 

Recalling that Q% (t) > Q^(t), Qk(t) > Q 3 (t) by construction, and using the uniform 
boundedness of Qjf and following from the uniform decay rate and the arguments 
above, we obtain for some M 

MW-oimpgir < 2M<mt)-o7{t))p$*,p8*>* 

||(Q S (*) - Qk{t))x\\ 2 < 2 M < (Q k {t) -Q s (t))x,x > H . 
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Again using (2.7), we have 


< (Qk (0 — Q ^ (O)-Ftf x > Pff x >H 

< < (Q^(h)-Q^(t t ))S N (t t ,t)P!jx,S N (t t ,t)P!jx >„, 

and 

< (Qk(t) ~ Q s {t))x,x > H << (Qk{tk) ~ Q s {t k ))S(t k ,t)x,S(t k ,t)x > H . 
Combining the above inequalities, we obtain 

IlCf (t)PS* - Q,(t)x \ | 2 < || Qi(t)P$x - Q k (t)x || 2 

+4j\/A/(||S A '(( l .«)i|r+||S(( l ,0if)- 

Let A; be large enough so that ||5(ffc, f)x|| < e/(12 MM), ||5' /v (tjt, t)Pffx\\ < t/(12M M), 
for all N. Then let N be large enough to obtain (t )P h x — Qit(f )x|| 2 — e /3. From 
the previous estimates we thus find \\Q^(t)Pff x — Q s {t)x\\ 2 < e which yields the 
desired results. 

We note that if Hypothesis 3.4 holds and the uniform decay assumption in 
Theorem 3.1 is replaced by the following: there exists t such that for any t > t and 
for any s 6 [to , */], \\S{t + s,s)x|| < e, 115^(1 + s,s)Pffx\\ < e, for all N, then the 
convergence of Theorem 3.1 is uniform in the bounded interval [f 0 , f/]. 

Theorem 3.1 is not very useful in practice, since the uniform decay assump- 
tion is difficult to verify directly. However it does provide some insight and sug- 
gests more realistic conditions that might be verifiable. Recalling the definition of 
W-detectability and our discussions following it, we are prompted to formulate the 
following assumptions. 

Hypothesis 3.6 (Uniform W-detectability) The original system is detectable 
and there exist constants Mk, Mkv and 0 > 0 such that for each N , there exists a 
uniformly bounded operator valued function K N {-) : H N h- >• H N , with \\K < 
Mk, for t € [to, oo). IfT% N is the evolution operator corresponding to the perturba- 
tion ofT N by K N VW", then \\T” N {t, s)\\ l{hN) < M Kv e~^- s \ for {t,s) € A^io). 
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If the Hypothesis 3.6 holds, then is the unique uniformly bounded solu- 
tion of the Riccati integral equation on the infinite time interval for H N . Under the 
uniform W-stabilizability Hypothesis 3.3, we have Qmin(t) ^ M-7, for all t G [to, oo). 
Furthermore, by application of Theorem 2.2, there exist constants M s ,a > 0 inde- 
pendent of N such that the evolution operator S N defined via Q^ in satisfies 

\\S N (t,s)\\ < M 3 e- a «~ 3 \ for (t,s) € A 0o (f 0 ). 

Thus, by Theorem 3.1, Qminif ) converges to Q m i n {t) as N — ► oo. We summarize the 
results in a major convergence theorem. 

Theorem 3.2 Assume that our system and its approximate systems satisfy Hypothe- 
ses 3.1 and 3.2 and the uniform W-stabilizability and uniform W-detectability condi- 
tions of Hypothesis 3.3, 3.6. Then the unique uniformly bounded solution Q N of the 
Riccati integral equation on [to, oo) in H N converges strongly to the unique uniformly 
bounded solution Q of the Riccati integral equation on the infinite time interval in H . 
Furthermore, if Hypothesis 3.2 is replaced by Hypothesis 3.4, then this convergence is 
uniform in t for t in any bounded interval. 

We note that in the case of a periodic system, the uniform convergence in one 
period implies that Q N converges to Q uniformly in the entire interval [to, oo). We fur- 
ther remark that the convergence of the Riccati operator guaranteed by Theorem 3.2 
is sufficient (using standard arguments, see [G], [BK 1] ) to guarantee convergence of 
the optimal approximate feedback system trajectories S N (t,s)Pf[x and optimal ap- 
proximate controls u N to the optimal system trajectories S(t , s)rr and optimal controls 
u (see Theorem 3.1 of [BI< 1] ) . Moreover, one also obtains convergence of the system 
generated by using the approximate feedback gains with the original infinite dimen- 
sional control system (a feature that is of great practical importance), e.g., see the 
related remarks in section 4 of [BK1]. 

The hypotheses of Theorem 3.2 are much more readily verified than others 
guaranteeing convergence that can be found in the literature (e.g., see [G], Theorem 
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5.3, where one is required to show that the approximate systems are uniformly sta- 
bilized by the feedback with a uniformly bounded sequence of approximate Riccati 
operators). As we shall see in the later sections, there are two distinct approaches 
that lead to rather easy use of our Theorem 3.2 in the event one is dealing with 
parabolic evolution systems. 

4 Parabolic evolution equations: control and approxima- 
tion 

In this section, we formulate the linear quadratic regulator problem for an abstract 
parabolic control system. We focus our attention on systems associated with a time 
dependent sesquilinear form. First we review the theory of parabolic evolution equa- 
tions (relying heavily on [T]) and extend some related results in a form applicable to 
control problems. Then a control system is defined for which general assumptions of 
stabilizability and detectability are made. A framework for approximation schemes 
is presented and conditions for convergence of the operators involved are discussed 
under assumptions of uniform stabilizability and uniform detectability for the ap- 
proximate systems. Our discussions here are in the spirit of the approaches taken in 

[BKl], |BI1], [BI2]. 

Let //, V be two complex separable Hilbert spaces with < v>h, < v >v as 
inner products and || • ||//, || • ||v as norms respectively. Let V * be the dual space of V 
with < •, • >v*,v denoting the duality pairing. The space V is assumed to be densely 
and continuously embedded in H, and thus there exists a constant c such that for all 
v> e f, iiV’ii// < cHt/’Hv'. Since for each element of H , we can define a bounded linear 
functional on F by < <p, ip >h , for rp € F, we have the usual embedding relationship 
V CH CV\ 

For each t in the interval [£o, oo), consider a sesquilinear form a{t\ •, •) defined 
on F x F. We assume throughout that a has the following properties: 

Hypothesis 4.1 (F-Continuity) For each bounded interval [t 0 , *i], there exists a con- 
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stant C\ such that 

(4.1) |or(t;y>,^)| < C|||v»||vWv, fort e[t 0 ,ti]^^ eV. 

Hypothesis 4.2 (V-Ellipticity) For each bounded interval [t 0 , t x ], there exist con- 
stants c 2 > 0, m such that 

(4.2) Re a(t;<p,<p) > c 2 \\<p\\v ~ ™\\<p\\h, for t € [t 0 , t x ], ip G V. 

Under the above assumptions, we have a well known ([FM], [K], [T], [S]) result: 
For each t (E [to, U], there exists a unique closed operator A(t) : V •-> V* such that 

(4.3) a(t-,<p,ip) = - < A(t)(p,xp >v\v, for xp 6 V. 

Furthermore, if A(t) is defined using the same method with a sesquilinear form a* 
defined by <r*(t; ip, ip) — <r(t; ip , <p), then A(t) is identical to the adjoint operator A*(t ) 
of A(t). Both operators A(t), A(t)* are infinitesimal generators of analytic semigroups 
in V* and an abstract parabolic evolution equation can be defined by 

^x(t) = A(t)x(t), x{t 0 ) = x 0 £ V. 

at 

In order to insure the existence of an evolution operator for this equation, we must 
make additional assumptions on the continuity of a with respect to t. 

Hypothesis 4.3 (Smoothness in t ) For each bounded interval [to, ti], there exist 
constants K and a, 0 < a < 1, such that for all t,s G [to, fi], and for all ip,tp € V , 
we have 

W(t]<f,ip) - cr(s] ip)\ < K\t - -ribllvIl^Hv. 

Under the above assumptions, there exists an evolution operator associated 
with the above evolution equation. The following theorem summarizes the properties 
of this evolution operator. 

Theorem 4.1 ([T], pp.127, pp. 145-155) Let the Hypotheses 4-1, f.2, 4-3 hold. Then 
there exists a unique evolution operator in V * satisfying the following condi- 

tions: 
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(i) For any t 0 < s < t < ti, the range TZ(T(t,s)) of operator T(t,s ) is a subset of 

V. 

(ii) The operator T(t,s)A(s) has a unique bounded extension in L(V*), for all t 0 < 
s < t < t\; therefore, we can and will use the same expression for the extension. 

(iii) For each <p G V*, the V* -valued function T(t,s)<p is continuously differentiable 

in t for t G (-Mi], and continuously differentiable in s for s G Further- 

more, for ip G V* 

^ T{t,s)p = -T{t,s)A(s)p. 

(iv) The restriction ofT(t,$) on H is strongly continuous in the H norm. For all 
xo G H, the function x(t) = T(t,s)x o is in L 2 ([s, ti] ; V) and the derivative 
x(t) = A(t)T(t,s)x 0 is in L 2 ([s,ti];V*). Furthermore, there exist constants 
Ci,C 2 , depending only on ci,c 2 ,m, K and a such that 

(4.4) ||f’(<,5)xo||v < Cfft - 5)- 1/2 ||x 0 ||//. 

(4-5) ||r(-,a)ar 0 ||La([. l t 1 ];V) < C 2 \\x 0 \\ H - 

All the statements in the above theorem can be found in [T]. However they are 
organized into several sections with somewhat different notation; we therefore give a 
brief argument which collects the results from the book. 

Proof: Existence. Taking X = V * , we let Aft) be defined as in (4.3). As indicated 
in [T], (pp. 144), using Theorem 5.2.1 of [T], we find there exists an evolution operator 
T on V*. The range TZ(T(t,s)) is a subset of V(A(t)) = V for all t 0 < s < t < t\. 
For any <p G V*, T(t,s)p is continuously differentiable in t for t G (^, tj]. Now let the 
sesquilinear form a* be defined by 
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Let A*(t ) be the linear operator defined via cr* ; then A*(t) is the adjoint operator 
of A(f). As indicated by the remarks following Lemma 5.4.6 of [T], we can use the 
results of section 5.2 (with S and A* replacing U and A of [T] ) to construct an 
operator-valued function S(t,s), such that for all t 0 < s < t < A*(s)S(t,s ) is a 
bounded operator in V *, and for any </> G V*, S(t,s)ip is continuously differentiable 
in s for s G [to, t). Furthermore, for <p G V* 

-^-S(t,s)<p = -A*(s)S(t,s)y>. 

In fact, S(t,s) can be constructed as follows. Let exp{tA*(s)} be the semi-group 
generated by A*(s) and we define 

S(t,s ) = exp{(f — s)A*(f)} + W(t,s), 

W(t,s) = exp {(r — s)A*(r)}f?(f,r)(fr, 

where the function R can be computed by iterative methods using 

R(t,s) - R 1 (rj,s)R(t,ri)dr) = Ri{t,s), 

with Ri(t,s) — ( A*(t ) — A*(s))exp{(t — s)A*(f)}. Then following the same type of 
arguments as in [T] (pp. 149), we can conclude that S(t,s) = T*(t,s). Therefore 
T(t,s)A(s ) has a unique bounded extension in V*. For all <p G V*, T(t,s)(p is strongly 
differentiable in s for s G [to, 0- 

Finally statement (iv) of the above theorem can be found in the sections 5.4 
and 5.5 of [T]. We note that in these sections of [T], the space X plays the role of our 
space H. Let T(t,s ) be the restriction of T(t,s ) to H] by Theorem 5.4.1 of [T], T(t,s ) 
is strongly continuous in the H norm. Furthermore the estimate (4.4) holds. For any 
x 0 G H, let x(t ) = T(t,s) xq. By Lemma 5.5.2 and Proposition 5.5.1 of [T] (pp. 152 
- 153) with / = 0, the function x(-) is in T 2 ([s,tj]; V) and x(-) is in L 2 ([s, <i]; V'*). 
In addition, the estimate (4.5) holds. We note that the constants C\,C 2 depend only 
on the constants ci,C 2 ,m,/f and a. 
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Uniqueness. By Theorem 5.2.3 of [T] (pp. 128), we can conclude that the 
evolution operator satisfying the conditions (i) - (iv) must be unique. 

We remark that the same theorem holds if we use the sesquilinear form a*; 
therefore, properties (i) - (iv) hold for the adjoint evolution operator T*(-,-). As 
consequence of (iv), the restriction T of T to H is an evolution operator in H as 
defined in the section 2. We wish to take H as our state space; therefore we use 
primarily the evolution operator T in this paper. The operator T is used in the 
remainder of the current section in several proofs of uniqueness theorems. The only 
precaution one must take is that T{t,s)<p is continuously differentiable with respect 
to t in the V * sense and the derivative of T(t,s)<p is an element of V*. In particular, 
for each xf E V, < T(t,s)<p,xf >// is differentiable with respect to t , and 

< T(t,s)<p,xp >//=< A{t)T(t,s)<p,xp >v*,v= -a{t-,T(t,s)tp,xl>). 

The conclusions of this theorem are very useful in defining our control system; 
however, the conditions of Hypothesis 4.3 are too restrictive for our use, since we may 
need to perturb the equation with nonsmooth but bounded (feedback) terms. We can 
show that if cr is perturbed with a sesquilinear form that is uniformly bounded in H, 
then there exists an associated evolution operator T k which preserves most of the 
desirable properties of the evolution operator T. In fact, let K(-) : [t 0 , oo) ► L(H ) 
be a uniformly bounded operator valued measurable function. We can then define a 
sesquilinear form ok in V x V as 

cri<(t;v,xl>) = cr(t]<p,xp)- < K(t)<p,if> > H , V, tp G V. 

It is easy to see that for each bounded interval Hypotheses 4.1 and 4.2 hold. 

Therefore, we can find an operator A^-(t) defined on V such that (4.3) holds for 
cr K and A^-(t). Furthermore, we have by the definitions of A(t) and Aj<(t) that 
Axil) 1 ? — A(t)tp -f K(t )<p for (p € V and we may establish the following result. 

Theorem 4.2 Consider a sesquilinear form a satisfying Hypotheses 4.1, 4-2, 4-3 and 
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let K{-), be defined as above. Then there exists a unique evolution operator Tk(-, •) 
in H for which the following properties hold: 

(i) The range 7£(Tr-(<,.s)) of the operator Tx(t,s) is a subset ofV, for all to < s < 
t < t x . 

(ii) For (p G H, the function TK(t,s)<p is differentiable with respect to t in the V* 
sense, and 

—Tx(t,s)ip = A K (t)T K (t,s)<p. 

(iii) For all xo G H , the function x(t) = Tfc(t,s)x o is in L 2 ([s,^i]; V ) and its deriva- 
tive x(t) is in L 2 ([s,ti]]V*). Furthermore, there exists constants C X ,C 2 depend- 
ing only on Ci,C 2 ,m, K and a such that 

\\T K (t,s)z 0 \\v < C x (t - s)- 1 ' 2 \\x 0 \\h, 

||r/r(-,5)x 0 ||L»([.,tiJ;V) < C 2 \\X 0 \\«. 

Proof: Existence. Let Tk be the unique evolution operator in H corresponding to 
the perturbation of T by K. From the results on perturbations given in section 2, we 
have that Tk satisfies for all <p G H 

(4.6) T K (t,s)ip = T(t,s)(p+ [ T(t,Tj)K(rj)T K (T],s)(pdT ?, 

J S 

T K (t,s)<p = T(t,s)<p + jf T K (t,ri)K(7j)T(T],s)^dg. 

Since the function &)<? IS uniformly bounded in H norm by some constant 

C, using the estimate (4.4), we can find a constant C such that for rj £ [s,f] 

Therefore the integral term in equation (4.6) converges in the V sense and hence, 
Tx(t,s)tp G V for ip G H. 

Consider x 0 G H, we define x{t ) = T^(t,t 0 )x 0 , and f(t) = K(t)x(t). From 
equation (4.6), the function x(t) can be written as 

x(t) = T(t,t 0 )x o+ f T(t, rj)f (g)di]. 

Jt o 
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By the strong continuity of Tk and uniform boundedness of K , it is obvious that 
/(•) E L 2 ([f 0 ,<i];tf) and hence /(•) (E T 2 ([f 0 i <i] ; V*). By Theorem 5.5.1 of [T], x(-) 
is in L 2 ([t 0 ,ti]-,V), is differentiable with i(-) in L 2 ([t 0 ,ti]-,V*) and satisfies x(t ) = 
Ax(t)x(t). Using the equality (4.6) and the boundedness of the perturbation, by 
modifying the constants Ci, C? in (4.4), (4.5), we can easily obtain 

||^(^o) < < 5 lk 2 ([to,<i];V) < C 2 \W\\h, 


for all <p €. H. 

Uniqueness. Let Tk satisfy the conclusions (i) - (ii) of Theorem 4.2. For all 
£ H, consider TK(t,s)<p as a V * valued function. Then we have 

^ T(t,rj)f K (ri,s)<p = -f(t,T])A(Tl)f K (T},s)ip 

+T(t, r,)(A(r,) + K(r,))f K (ri,s)<p 
= f(t,T])K(T))f K {Tl,s)ip. 


Integrating both sides of the above equation from s to t, we obtain 

Th(t,s)<P = T(t,s)<p + £ f(t,Tj)K(j))f K (r},s)(pdri. 

Since T is the restriction of T to H, Tk is a solution of (4.6). By our uniqueness 
results of section 2 for perturbed evolution operators, we have Tk = Tk . Hence the 
unique solution of (4.6) is the unique evolution operator Tk generated in the theorem. 

Now consider a function /(•) € L 2 ([t 0 ,ti\; H). We can then define 

z(t) = T(t,t 0 )z 0 + j T(t,T])f(Tj)dr]. 

Jt o 

The function z{-) is the unique solution of the following initial value problem: 

(4.7) z(t) = T(t,s)z(s) + J' T(t,r))f(rj)dr), t 0 <s<t<U. 

z(t Q ) = z 0 . 
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Henceforth, we consider (4.7) as the definition of our basic evolution system. The 
function z corresponds to the solution of a weaker formulation of the evolution equa- 
tion. 

Lemma 4.1 ([T], Theorem 5.5.1) The function z(-) given by (4-7) is the unique 
function in L 2 ([to, <i]; V) with derivative i(-) in L 2 ([t 0 , H]; V*) for which the following 
equation holds for if € V 

(4.8) < z(t) -z(t 0 ),il> >h = [ {-°(r, z (v ),V0+ < /M’V’ >//}<*?, 

Jto 

z(to) = 2 ' 0 - 

Lemma 4.2 ([T], Lemma 5.5.1) For any two functions 2 (-),u;(-) in L 2 {[t 0 , <i]; V) 
with derivatives z,w in L 2 ([to, $i]; V*), the following equality holds: 

< z(t),w(t) > H = <z(s),w(s)> H 

+ I {< z{rj), w(g) >v\v +< w(T]),z(y) >v\v}dq, 

J s 

for all t 0 < s < t < t\. 

As a consequence, for any a;o € H, let x(t) = T(f,<o) ar o, then 

(4.9) lk(0llff = IMItf - 2 / Re a(y;x(T]),x(Tj))dq. 

Jto 

We note that if Hypotheses 4.1 - 4.3 hold, then for each bounded interval 
[f 0 , <i], we can define T(t,s) uniquely, therefore T(t,s) is also uniquely defined for all 
t 0 < s < t < oo. The equality (4.9) suggests a sufficient condition for the stability of 
T. 

Hypothesis 4.4 There exists a constant k > 0 such that 

Re <r(t;tp,ip) > k\\tp\\ 2 H , fort 0 <t< oo, ip € V. 

Theorem 4.3 Under Hypothesis 4 - 4 > T is uniformly exponentially stable. 
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Proof: For any xo G H , let x(t) = T(t,s)x o- Then by (4.9), we have 

IMOIIh ^ IMI»-2jf *||®(i?)|| 2 <fy, 

for all t 0 < s < t < oo. This implies 

l|r(M)lk(«) < 1, J‘ m-M)*oll« < 

Therefore, by Lemma 2.1, T is uniformly exponentially stable. Note, moreover that 
by Lemma 2.1, under Hypothesis 4.4, we can find M,a > 0 depending only on k such 
that 

linM)ll«H> < Me'" 

We can now use these considerations to define an evolution equation con- 
trol system of the form (4.7) via a sesquilinear form. The space H will serve as 
our state space, with subspace V and the sesquilinear form a defined as above and 
Hypotheses 4. 1-4.3 holding. Let the control space U be a Hilbert space, and let 
B(-) : [to, oo) t-* L(U, H ) be a measurable operator-valued function. We assume that 
there exists a constant Mb such that 

II B(t) || l(u> h) < M b , for i G [t 0 , oo). 

For any control rf(-) : [to, oo) i— »• U, belonging to L 2 ([t 0 , oo); U), the corresponding 
trajectories satisfy for tp G V 

(4.10) < z(t) - z(s),ip > H = - £ {(t(t);z(t)),iI>)- < B(t})u(tj),^ > H }dr), 

for all (t,s) G Aco(to)- Let be the evolution operator defined via a. By 

Lemma 4.1, an equivalent form of (4.10) is given by 

(4.11) z(t) = T(t,s)z(s) + £ T(t,r))B(Ti)u(ri)dr}, for (t,s) G A oc ,(t 0 ). 

Let zo G H be the initial state of the system at to and let the cost for control u(-) be 

given by 

roo 

(4.12) Joo{w, zo,to) = / < W(t)z(t), z(t) > H + < R(t)u(t),u(t) >u dt, 

Jto 
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where W(-) : [t 0 , oo) i— > L(H), R(-) : [t 0 , oo) t— ► L(U). The operators W(t),R(t) are 
assumed to be selfadjoint nonnegative definite operators, uniformly bounded in the 
entire interval [t 0 , oo). Furthermore, there exists a constant r > 0 such that 

<R(t)v,v>u > r||u||y, for t > t 0 , v £ U. 

Recalling the discussions of section 2, we note that the standing assumptions 
of that section hold. Therefore, for a given nonnegative definite self-adjoint operator 
G on H, the Riccati integral equation in each finite time interval [to, <*], 

Qk(s)x = T*(t k ,s)GT(t k ,s)x 

+ /“ T- [v W(v) - QtMBMR-'MB-MQ.M] T(v,s)xd v 
has a unique self-adjoint solution Q k . 

For the control problem in the infinite time interval [to, oo), we need stabiliz- 
ability and detectability conditions to assure existence and uniqueness of a uniformly 
bounded solution of the Riccati integral equation. 

Hypothesis 4.5 (Detectability) There exists a uniformly bounded operator valued 
function ^(-) : [t 0 , oo) L(H), such that if we denote by 5$ the evolution operator 
corresponding to the perturbation of T by 'f r (-)fF 1 / 2 (-), the following estimate holds 
for x 6 H : 

\\S*(t,s)x\\ H < Me-^\\x\\ H , 

for some constants M,u> > 0. 

Hypothesis 4.6 (Stabilizability) There exists a uniformly bounded operator valued 
function K(-) : [t 0 , oo) L(H , U), such that if we denote by Sk the evolution operator 

corresponding to the perturbation of T by B(-)K(-), the following estimate holds for 
x G H: 

\\S K (t,s)x\\ H < Me-^\\x\\ H , 

for some constants M,u > 0. 
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We remark that Hypothesis 4.6 is stronger than “W-Stabilizability”; however 
under the Hypothesis 4.5 by Theorem 2.2, these two types of stabilizability assump- 
tion are equivalent. 

To this point we have defined a control system using an abstract parabolic 
evolution equation that fits into the general framework of section 3. Under Hypothe- 
ses 4.5 and 4.6, we may apply the theory of the previous sections to establish the 
following results for our control problem: 

(i) The Riccati integral equation in the infinite time interval [t 0 , oo) has a unique 
uniformly bounded solution Q(-). 

(ii) Let Sq be the evolution operator corresponding to the perturbation of T(-, •) by 
-BR~ 1 B*Q(-). For each initial state zq , the unique optimal trajectory is given 
by S Q {t,to)z 0 . 

We turn next to give results for finite dimensional approximations of our con- 
trol system. As in section 3, let be a sequence of finite dimensional sub- 

spaces of V C H. Let P$ be the orthogonal projection operator from H onto H N . 
Since H N is an approximation of i/, we assume that for every v? G H, WP^^ — ^Wh — ► 
0, as N —> oo. In addition, we require that H N is an approximation of V as well, 
so that for all y € V, || P$<p — <^||y -+ 0, as N — ► oo. We note that in fact this 
latter convergence implies the convergence in H for (p € H since V is continuously 
and densely embedded in H. 

Let {U N }'^-. 1 be a sequence of finite dimensional subspaces of U. Let Py be 
the orthogonal projection operator from U onto U N . We assume U N approximates 
U in the following sense: for v 6 U, \\PyV — v\\y — ► 0, as — > oo. 

For each N , we define a sesquilinear form cr N as the restriction of cr to H N x H N 
and define a linear operator A N (t) : H N i— ► H N by 

- < A N (t)<f N ,iJ> N > H = a N (t-,<p N ^ N ), for € H N . 

By continuity of a with respect to t, the operator valued function A N (t ) is continuous 
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in time. As a consequence, there exists a unique differentiable evolution operator 
r"( v ) in H N generated by A N (t)\ that is for <p N € H N we have 

jT K (t,s)<p K = A N (t)T N (t,s) V N . 

Note immediately that the a N ' s satisfy the Hypotheses 4. 1-4.3 with the same con- 
stants cj, c 2 , m, a and K, therefore for each fixed interval [t 0 , U], there exist constants 
C\, C 2 independent of N such that for all <p N € H N , 

(4.13) l|T''(<,sV' V ||K<C,((- S )-i|| ¥ .' v || wl 

(4.14) < C 2 \\ V n \\ h . 

The approximation properties of the evolution operator T N are summarized by the 
the following convergence theorem. 

Theorem 4.4 Let Hypotheses 4-1- 4-3 hold and let T(-,-) and T N (-,-) be defined as 
above where \\Pfj<p — y>||v — * * 0 as N — ► 00 for p> £ V; then the following properties 
hold: 

(i) There exist constants Mr and ui such that for all N , 

r' v (M)IU ( „» ) < M T e“ <•->, ||r((,s)||t( H ) < M T e*'-\ („<»<(< oo. 

(ii) For any finite interval [a, b] C [to, 00 ) and any ip 6 H , we have 

\\T N (t,s)Pff(p — — ► 0, asN — ♦ 00 , a < s < t < b. 

Furthermore the convergence is uniform for all a < s < t < b. 

Proof: (i) By Lemma 4.2, and (4.9) for every <p £ H N we have 

ll^*, s )t\\h = IMItf ~ 2 j s Re <r N {w tN (v y 8)<p, T N (r), s)<p)dy. 

Under Hypothesis 4.2, 

l|r"(MMIjr < Ml + 2j‘m\\T N ( n ,sM\ld v . 
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Using Gronwall’s inequality, we obtain 


< IMIte 21 ” 1 ''-'’- 

Noting that the same estimates hold for T(t,s)<p, we obtain (i). 

(ii) Let <p E H, define w{t) — T(t,s)ip and w N (t) = T N (t,s)Pjj <p and let 
z N (t) = w(t) — w N (t). We note that z N (t ) is not an element of H N , in fact 

(4.15) z N (t ) — P%z N (t) = w(t) — Pgw(t). 

Since w(t) is differentiable in the V* sense, w N (t) is differentiable, and both func- 
tions are in L 2 ([a,b]\V) with derivatives in L 2 ([a,b];V*). By (4.9), Lemma 4.2 and 
definitions of the operators A(t), A N (t), we obtain 

= ||A0)|| 2 h + 2 / Re < MvWl) - A n (t])w n (t]), z n (tj) >v,v dr) 

= \\z N (°)\\h - 2 1 Re {°(v™(v),z N {v))-v N (mw N (Tl),P%z N {Ti))}dTi 
~ 2 J s < A N (T])w N (ri),z N {Ti) - P%z n (t] ) >v,v dr). 

Since the duality pairing reduces to the H - inner product on H x H, we have 

<A N (r,)w N (r,),z N (r,)-P>lz N W >y.y 

=< A n Ww n W,z n W - >„ . 

Moreover, P$ is the orthogonal projection operator and hence the last term in the 
above equation equals to zero. Using the definition of cr N , the sesquilinearity of a 
and (4.15), we find 

\\z N {t)\\ 2 H = \\v-PHV\\ 2 H~ 2 J a te<r(riiw(ii) t z N (r))-Pgz N (Ti))dri 

-2 1 Re o{r)\ w(rj), P%z N (rj))dr) 

+2 L Re aN (VA t}N {v),PH zN { T l))dv 
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= \W - P H <p\\h - 2 J s Re a (v, Hv), w(v) - P H w{rj))dr] 

-2 jf Re a( rj] w(r)) - P% w(rj), P% z N (T)))dj] 

~ 2 I s Re 

Since Pjj is the orthogonal projection, we find < PjjZ N ,w — Pfj w >H = 0, so that 
from (4.15) 

\\pg^m 2 n = \\z N m 2 H - iko - pg*(t)\\H- 

Combining this with the previous equation, we have 

|K^(<)II« = Q N (t,s)-2j\e <, N ( r ,P%z N ( v ),P$z N (r,))d^ 
where 0 N (t,s) is given by 

e N (M) = \\<p - p%ip\\ 2 H - IMO - PhH^Wh 

rt 

-2 J Re a(r];w(r}),w(r}) - P%w(r]))dr) 

-2 jf Re a(i i; io(r j) - P% w(rj), P% z N (r]))dT). 

Using the U-ellipticity of cr^, we find 

(4.16) l|C£' 2 ' v (i)llt < \&"(t,s)\ + 2|m| j‘ \\PHz N M\\Ul- 

To use Gronwall’s inequality to conclude convergence of Pfj z N , it suffices to show 
that |0^(<,s)| goes to zero uniformly for all a < s < t < b. By the continuity and 
uniform boundedness of T, the term \\T(t,s)(p — PffT(t,s)ip\\ 2 H goes to zero uniformly 
for all a < s < t < b. Using the U-continuity of <j, the two integrals in 0 can be 
bounded by 

2c i / {\\ p h w (v) - w(v)\\v\\Phz n (v)\\v + IM »?) IMI p h w (v) - «»(*/) l | v } d v 

v 5 

I 

< 2c, [|M-)|| t > (MllV0 + [jf IHf) - PXwmUn\- 
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Using the inequalities (4.5) and (4.14), we observe that the functions w,w N are in 
a bounded subset of L 2 ([a,b\-,V). By dominated convergence arguments, the above 
integral converges to zero. Furthermore by taking t = 6, s = a, we obtain that this 
convergence is uniform for all a < s < t < b. Therefore |0 /v (t,s)| converges to zero 
uniformly for all a < s < t < b. Finally, from (4.15), we have 

||T(i, s v- r w (M)/#dli = 11^(011?, + MO - C?«>(0 lit. 

and the uniform convergence of Pjf z N (t) in t implies T N (t,s)P N ip converges to 
T(t,s)<f uniformly for all a < s < t < b. 

Since we can define operators A*(t),A“ N ,T*(t,s) and T*(t,s) by using the 
sesquilinear form cr* as we indicated after (4.3), the convergence of T* N (t,s) to T*(t, s) 
can be shown using the same arguments as in the proof of the above theorem. 

Having defined our approximate (uncontrolled) system and established the 
convergence of Theorem 4.4, we return to the control problem for (4.10)-(4.12). Ap- 
proximations of functions B, W, R are defined as follows: 

B n {-) : [t 0 , oo) i — ► L(U N , H N ), B N (t)v N = P$B(t)v N , v N G U N ; 

VF A '(-) : [t 0 , oo) L(H n ), W N (t)<p N = P%W(t)v N , <f N G H N ; 

R n (-) : [t 0 , oo) i — ► L(U"), R N {t)v N = P$R(t)v N , v N G U N . 

Let G be the nonnegative selfadjoint operator in the finite interval cost functional 
associated in the usual manner with (4.12) for our control system in H. Let G N = 
P$G and z* = P$z 0 . 

In each subspace H N , a finite dimensional control system is thus defined by 

(4.17) z N (t ) = T N (t,s)z N (s) + J‘ T N (t, V )B N ( V )u N ( V )d V , 

with u N (‘) G L 2 ([to, oo); U N ) and ^^(to) = Zq . The cost functionals for the associated 
finite time interval problems are given by 

(4.18) J N (xi N -,z”,t 0 ,t k ) =< G N z N (t k ),z N (t k ) >„ 

+ /“ {< W N (t)z"(t),z N (t) > H + < R K (t)u N (t),u N (t) > v \it , 

J to ^ 
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while the cost functional for the infinite time interval problem is given by 


(4.19) jg(u N -,zg-,t 0 ) = 

Jf{< W N (l)z"(t),z"(l ) >„ + < R”(t)u N (t),u N (t) > v }dt. 

To obtain the uniform convergence of the operator valued functions, we make 
additional assumptions on the continuity of B, IT, R. 

Hypothesis 4.7 (Parameter smoothness) The operator valued functions B , W , R are 
piecewise strongly continuous functions on [t 0 , oo). 

Lemma 4.3 Under Hypothesis J t . 7, the following convergence is uniform in t for t 
in any bounded interval: 

II B N (t)Pffv - B(t)v\\„ — *■ 0, veU ; 

\\B* N (t)Pff(p — B*(t)tp\\u — » 0, 

\\W N {t)Pffy-W{t)y>\\ H ->Q, € H; 

|| R N {t)P$v - R(t)v\\u -* 0, veU, 

as N — > oo. The operator G N Pff also converges strongly to G as N — * oo. 

Proof: We only prove the uniform convergence of B N , the remainder of the arguments 
being similar. For simplicity, we without loss of generality assume that the function 
t i — ► B(t ) is strongly continuous. For a given v € U, the pointwise convergence of 
the functions B N (t)Pffv to B{t)v is given by our assumptions on the approximation 
properties of the spaces H N , U N \ To conclude uniform convergence in it is enough to 
show that the functions B N (t)PjjV are equi-continuous. That is, for any e > 0, there 
exists 6 > 0 such that for all N, if \t — s| < 6, we have \\B N (t)Pff v — B n (s)P/j v\\h < 
By definition of B N , we have 

\\B N (t)Pgv - B N {s)Pffv ||„ < || B(t)Pgv - B{s)Pgv ||„ 

< || B(t)(Pgv - r)||„ + 115(1)0 - B(i)»||» 
+l\B(s)(Pgv - o)||„ 

< 2M B \\Pgv - viler + \\B(t)v - B(s)v\\„. 
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By continuity of B , we conclude that B N (t)PyV are equi-continuous functions of t in 
a bounded interval. Hence the convergence is uniform in any bounded interval. 

It is easy to verify that B N , W N , R N are uniformly bounded and W N , R N are 
nonnegative self-adjoint operators. In addition there exists a constant r > 0 such 
that for all N, 

< R N (t)v N ,v N >u > r||u N ||^, for t G [to, °°), v N 6 U N . 

Consider any finite time interval [to, tk], and let Qk,Q k be the unique self- 
adjoint solutions of the Riccati integral equations in H and H N associated with the 
control systems (4.11), (4.17) respectively. Then it follows from Theorem 4.4 and the 
discussions of section 3 (in particular see Hypotheses 3.1, 3.2 and the remarks just 
prior to Hypothesis 3.4) that for each k, Q k {t) converges to Q k (t ) strongly and the 
convergence is uniform in t for t € [to, t*]. 

We assumed above (Hypotheses 4.5, 4.6) that the control system (4.10)-(4.12) 
in H is detectable and stabilizable. Therefore there exists a unique uniformly bounded 
solution Q of the Riccati integral equation on the infinite time interval [t 0 , oo). In 
order to approximate Q by a uniformly bounded solution of the Riccati integral 
equation in H N , we have to show that the approximate control systems defined here 
are also detectable and stabilizable. More importantly, recalling the results (e.g., see 
Theorem 3.2) of section 3, we need uniform detectability and uniform stabilizability 
for the approximate systems (4.17), (4.19). 

Based on the stabilizability and detectability properties of the original system, 
for a given approximation scheme, we would like to show that the following conditions 
hold: 

Condition US (Uniform stabilizability) There exist constants Mk,M, u > 0 inde- 
pendent of N such that for each of the approximate systems, we can find a uniformly 
bounded operator valued function K N {-) ■ [<o, oo) e- ► L(H N ,U N ) such that 

||-ft' N (£)||L(tf N ,l/ w ) < Mk, 
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and ifTff is the evolution operator corresponding to the perturbation ofT N by B N K N (-), 
then 

\\TlUt,s)\\L(H’‘) < A/e-"**-', for (t,s) € A «,((„). 

Condition UD (Uniform detectability) There exist constants > 0 indepen- 

dent of N such that for each of the approximation systems, we can find a uniformly 
bounded operator valued function '5 N ( •) : [f 0 , oo) i-* L(H N ) such that 

^ My, 

and if T$ is the evolution operator corresponding to the perturbation of T N by 
then 

l|r»((,s)||t ( H») < Me"""-), for ( t,s ) 6 A 

We may summarize our findings as follows. 

Theorem 4.5 Under Hypotheses 4-5~4-7, the conditions \\Pfj<p — <^11^. —> 0 

for (p e V, \\P/jV — u||i/ — > 0 for v G U and the Conditions US and UD, there exists a 
unique uniformly bounded solution Q N of the Riccati integral equation on the infinite 
time interval [£ 0 , oo) for each approximate system in H N . Furthermore, the sequence 
Q N {t)Pfj converges strongly to Q(t) and the convergence is uniform in t fort in any 
bounded interval. 

These results follow from Theorem 3.2 and our discussions above. We have 
thus reduced our problem of ensuring convergence of the Riccati variables to one of 
guaranteeing uniform stabilizability and detectability of the approximate systems. In 
the following two sections, two different approaches to obtaining uniform detectability 
(Condition UD) and stabilizability (Condition US) are presented. 

5 Dissipativity and uniform stabilizability / detectability 

The original control system (4.10) defined in section 4 was assumed to be stabilizable 
and detectable (i.e., we assumed Hypotheses 4.5, 4.6 held). For a given evolution sys- 
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tem, often an easy way to ascertain stability is using the dissipativity of the system. 
In particular, if a system satisfies Hypothesis 4.4, by Theorem 4.3 the associated evo- 
lution operator is uniformly exponentially stable. This naturally suggests a sufficient 
condition for stabilizability of a control system. 

Hypothesis 5.1 There exists a uniformly bounded function K(-) : [< 0 , oo) >-> L(H , U) 
and a constant k > 0 such that for all <p G V, 

<r(t;<p,tp)+ < B(t)K(t)(p,tp> H > *|Mljjr> for t e [t 0 , oo). 

Lemma 5.1 Under Hypothesis 5.1, the control system defined by (4-10)-(4-12) is 
stabilizable. In fact if Tk is the evolution operator corresponding to the perturbation 
of the evolution operator T by —BK, we can find constants M, a > 0 such that 

||7k(M)||l ( *) < for (: t,s ) G A^U). 

As a consequence there exists a constant C such that for all x 0 £ H , we can find a 
control u(-) £ Z/ 2 ([s, oo); U) with a cost 

Joo (u-,x 0 ,s) < Cllxoll 2 , for s G [to, oo). 

Proof: Let K(-) be the operator valued function in Hypothesis 5.1. Define the 
perturbed sesquilinear form cr A - ( t ; p, xfi) = a (t-,<p,ij>)+ < B(t)K(t)<p,if > H . Then T K 
is associated with cr x as in Theorem 4.2 with —B(t)K(t) as the perturbation term. 
Under our assumptions, by Theorem 4.3, there exist M, a > 0 such that 

IIIMMlIkm < JVe-"«‘->, for (t,») € A 

For any x 0 G H, let v(t) = — K(t)Tfc(t,s)x 0 ; it is easy to see that the corre- 
sponding trajectory is x(t) = Tx(t,s)x 0 . 

By our standing assumptions, the operator valued functions W(-), R(-), B(-) 
are uniformly bounded in the entire interval [ t 0 , oo). Take 

C > (II W(*)||l(H) + \\K*m(t)K(t)\\ L(H) )M 2 /2a, for t > t 0 . 
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Then 


r co 

J<x,(v;x 0 ,s ) = J {< W(t)x(t),x(t) >h + < R{t)K(t)x(t),I<(t)x(t) >u} dt 

< cn*„r„. 

Similarly, a sufficient condition for detectability can be stated as following. 

Hypothesis 5.2 There exists a uniformly bounded operator valued function \P(«) : 
[t 0 , oo) L(H), and constant A > 0 such that 

¥>,¥>)+ < y(t)W*(t)<p,<p > H > A||v?||^, <p € V. 

Lemma 5.2 Under Hypothesis 5.2, the control system defined by ( 4 . 10)- ( 4-12 ) is 
detectable. 

In the remainder of the current section, we assume that Hypotheses 5.1, 5.2 
hold for our control system in H. The strict i/-dissipativity Hypotheses 5.1, 5.2 on 
the evolution systems are stronger than the usual stabilizability and detectability 
hypotheses; however, they are in general easy to verify for a wide class of problems. 
Moreover, the constants M and the decay rates a depend only on the values of k and 
A. Thus, this type of approach suggests that approximate systems which preserve the 
/f-dissipativity might be uniformly stabilizable and uniformly detectable. Pursuing 
this type of argument, we shall try to show that the following conditions are implied 
by Hypothesis 5.1, 5.2. (As in the discussions of section 4 surrounding (4.17)- (4.19), 
we assume that B N {t ) = Pff B{t) and W^(0 = PffW{t ) .) 

Condition 5.1 There exists a constant k > 0 such that for every N , there exists a 
uniformly bounded operator valued function K N (-) '■ [<o, °°) L(H N ,U N ) so that 

(5.1) <, N (f )V , N tV N )+ < B"(t)K N (t)Y S ,v N > H > t||?X 

holds for all <p N € H N . 
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Condition 5.2 There exists a constant A > 0 such that for every N, there exists a 
uniformly bounded operator valued function : [to, oo) •-* L(H N ) so that 

(5.2) a N (t;<p N ,V N )+ < >H > \\\<p N \\ 2 H 

holds for all tp N G H N . 

Note that if the original system satisfies Hypotheses 5.1 and 5.2, by the defi- 
nition of a N we have 

(5.3) ,‘p N )+ < Pff B(t)K(t)^ N >h > kW^Wl 

(5.4) a r, (f,V N ,V N )-KPSmwHt)f N , 9 N > H > A|| v > w ||?,. 

Let us compare inequality (5.3) to (5.1); if we could take K N (t) = K(t) in (5.1), 
then Condition 5.1 holds trivially. However, a careful examination of inequalities 
(5.3) and (5.4) reveals that they do not provide stabilizability and detectability of 
the approximate system. In the case of (5.3) vs. (5.1) we observe that the range of 
the operator K(-) hypothesized in Hypothesis 5.1 is not necessarily in U N and K(-) 
cannot be used as a stability operator for the control system in H N , U N as required 
in Condition 5.1. Comparing (5.4) and (5.2) and recalling that W N (t) = PjfW{t ), 
we see that the the choice = Pjf'ilPfj would suffice only in the case where 
Pfl(P%W(t)) 1/2 = 

Before we state additional conditions for the approximate systems, let us con- 
sider several interesting cases for which stabilizability and detectability are preserved. 

I) Dissipative systems. Suppose that Hypotheses 5.1, 5.2 hold for K{t) = 0 and 
'l'(t) = 0; then by definition of the sesquilinear form a N , the Conditions 5.1, 5.2 hold 
with K N {t ) = 0 and 'I ,N (t) = 0. This is the case when the homogeneous system is 
itself dissipative. 

II) Finite dimensional control. Suppose the control space U is finite dimen- 
sional. Taking U N = U, we can use K N (t) = K(t) and the approximate systems are 
uniformly stabilizable. 
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Ill) Special stability operators. Consider the inequalities in Hypotheses 5.1, 5.2 
again. We can assume that these inequalities hold where B(t)K(t) and 
are nonnegative definite self-adjoint operators. Suppose that there exist scalar func- 
tions K,(t) > 0 and fi(t) > 0 such that 

B(t)I<{t) < K (t)B(t)B*(t), <H{t)W 1/2 (t) < fi(t)W(t). 

Then taking K(t ) = K(t)B*(t) and \P(<) = ii(t)W x ! 2 (t), we find Hypotheses 5.1 and 
5.2 also hold. If we modify slightly the definition of the sesquilinear form <r N by 

+ <mv - p$)wHt)v N ^ 

the sesquilinear form a N satisfies the Conditions 5.1 and 5.2. Indeed we note that 
the perturbation terms satisfy 

<B(t)[I-Pg] K(t)<p N ,r N > 

= K(t) <[I- P$}B*(t)<p", [I - P#]B-(t)<p N >u > 0, 

< V{t)[I - P%]W ll \ N ,<p N > 

= f i(t) < [I - P»}W'' 2 (t)<p N , [I - P%)W'' 2 {t)y N > H > 0. 

Thus by taking I< N (t) = Pff I\(t), the Conditions 5.1, 5.2 hold for 

cr N • On the other hand, the perturbation terms go to zero as N goes to oo. Therefore 
if we use a N as the sesquilinear form for the approximate control system in H N , the 
corresponding evolution operator T N should also converge to T. 

The three cases above motivate us to consider the following modifications of 
the sesquilinear form in H N . Let the operator valued functions A'(-), 'P(-) be as in 
Hypotheses 5.1 and 5.2; define a N as: 

a N (t;^ N ^ N ) = a N (t;<i> N ,i> N )+<B(t)[I-Plf}K(t)<p N ,il> N > 

+ <*(*)[/ >, 
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for all <p N ,tf N £ H N . Let A N (t) : H N i— ► H N be defined by 

- < A N (t)<p N , r > H = a N (t ; <p N , if N ), <p N , £ H N . 

Let T N (-, •) be the evolution operator generated by A N (t). 

We can repeat the arguments in the proof of Theorem 4.4 using T N in place 
of T n . In the arguments, there is an extra term 

2jf Re < A n (t})w n (t}), z N (rj) > H dp 

on the right side of the inequalities, where A N (r/) 6 L(H N ) is given by 

A n (() = P$B(t)[I - Pff)K(t) + P8*W - P£}W''\t). 

There exists a constant C independent of N such that 

||A ;V ( 7 ?)||l(// jv ) 5: C, r] E [a, 6 ], 

and, furthermore, 

II A N {v)Ph v\\h -» 0, for (f £ H. 

Recalling z N = w — w N and using (4.15), we have w N — w — z N = PfjW — Pjj z N . 
We thus find (suppressing the argument tj throughout) 

Re < A V, z N > H = |< A N (Pgw - P%z N ), P"z N > H 

< \\* n pZ'°\\h\\p!I* n \\h + c\\pHz n \\ 2 h 

< \\\A N pxm+(c+ 1 -)\\pgz K ui. 

The integral (with respect to 77 ) of the first term in this last expression — ► 0 uniformly 
in t,s and can be added to the term O^^s) in (4.16), while the integral of the second 
term can be included with the integral term in the right side of (4.16). We thus can 
argue: 

Theorem 5.1 Under Hypotheses 5.1, 5.2, the conclusions (i), (ii) of Theorem f.f' 
hold for T N . 
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Now consider T N as the evolution operator for our approximate control systems 
in H n \ the convergence of the solutions of the Riccati integral equation in any finite 
time interval still holds. To generalize the arguments in the three special cases I), II), 
III) above, we make the following additional assumptions. 

Hypothesis 5.3 Consider /{(■), &(■) as in Hypotheses 5.1, 5.2 and assume there 
exist constants k < k, A < A and N such that for all N > N 

<(I - Pu)K(t)<fi N .(I - Py)B'(t)(fi N > > -MI'/'IIh. 

<(l-PS)wi(t) v N ,(I-P%W(t)v N > > -A||v>X, 

for all <p N G H n . 

Lemma 5.3 Under Hypotheses 5.1, 5.2 and 5.3, the approximate control systems are 
uniformly stabilizable and detectable. 

A A. A 

Proof: We assume without loss of generality that N = 1. Let k = k — k and 
A = A — A. By Hypothesis 5.3, k > 0 and A > 0. Take K N (t) = P/f K(t) and 
ty N (t) = Pffty^t)', then with this choice of K N ,ty N , the Conditions 5.1, 5.2 hold for 
a N . Let Tjf , T§ be evolution operators corresponding to the perturbations of T N by 
B N K N and ^ ;v (VR ;v ) 1 / 2 respectively. By Theorem 4.3, there exist constants M, a > 0 
depending on k , A only such that 

l|f£(M)ll«»», < Me - <•->, < Me-"-'. 

By the general framework of section 3, there exists a unique solution Q” of 
the Riccati equation on the infinite time interval for each control system in H N . 
The operator Q N (t)Pff converges strongly to the unique solution Q(t ) of the Riccati 
integral equation for the original system in H as jV — + oo. The convergence is uniform 
in t for t in any bounded interval. 

We summarize the results for our dissipativity approach to uniform stabiliz- 
ability and detectability in the following theorem. 
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Theorem 5.2 Consider the parabolic control system defined by (4. 10)- (4-12) under 
Hypothesis 4-l~4-3, 4 -5-4-7, and the corresponding approximate systems as defined via 
cr N , T N as in this section where P — * I strongly in V and Pff — ► I strongly in U . 
Under the “H-dissipativity” Hypotheses 5.1, 5.2 and the consistency Hypothesis 5.3, 
the following conclusions hold: 

(i) There exist unique uniformly bounded solutions Q N , Q of the Riccati integral 
equations in the infinite time interval [/ 0 , 00 ) for each of the approximate control 
systems and the original system, respectively. There exists a constant M such 
that for all t £ [to, 00 ) 

(ii) Let S N , S be the perturbed evolution operator corresponding to the perturbations 
of j>N t by — B N (R N )~ l B* N Q N and —BR~ 1 BQ respectively; then there exist 
constants M, a > 0 such that 

l|S W (<,»)lll.(H») < A/e-o' 1 -'; ||S((, S )|| lm < 

for all (t,s) G A co (< 0 ). 

(iii) As N —* oo, Q N ,S N converge to Q,S in the following sense: for all £ H 

\\Q N (t)PHv-Q(tM\H -><>■, 

l|s"(i, - S(<, <M|h -* o. 

The convergence is uniform in ( t,s ) in any bounded interval. 

The advantage of using the dissipativity approach outlined above is that the 
hypotheses are readily checked. The .//-dissipativity can sometimes be replaced by 
even weaker dissipativity conditions for which one can obtain an exponential decay 
rate (e.g., see [Ch], [La]). For parabolic systems with strict F-ellipticity, we can 
avoid use of this type of argument as we shall see in the next section. However these 
results might be useful for systems without strict V-ellipticity or possibly even some 
hyperbolic systems (e.g., see [BKS], [BKW]). 
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6 Periodic systems: compactness and uniform stabilizabil- 
ity / detectability 

One of the special features of parabolic evolution systems as defined in section 4 is 
that the evolution operator T(t,s ) is also a bounded linear operator from H to V. 
Since often the space V is compactly embedded in H , T(t,s ) is thus a compact oper- 
ator. Using this fact, we can show that the convergence of the sequence of operators 
T N (t,s ) to T(t,s) is in a stronger sense. In this section, by combining periodicity and 
compactness of the evolution operators T N and T , we can show that the approxima- 
tion schemes discussed in section 4 preserve detectability and stabilizability. 

Some ideas for the stability of periodic evolution operators can be found in 
[HI], [H2]. The use of compact embedding ideas for the proof of operator norm 
convergence can be found in [BI2] ( The authors gratefully acknowledge K. Ito for 
fruitful discussions regarding this approach). 

In this section we make a periodicity assumption for our control system: 
Hypothesis 6.1 There exists a constant 0 > 0 such that 

(i) The sesquilinear form <y is 0-periodic in time; 

(ii) The operator valued functions B, R , W are periodic in time with period 0. 

Lemma 6.1 Under the above assumption, the evolution operator T(-,-) of the corre- 
sponding homogeneous evolution equation is 0-periodic. 

Proof: For any s < t, and all € V, we have 

rt+e 

< T(t + 0,3 + >h - o{t ;T(t,s + $)v,%l>)dT 

Js+6 

= < >H -J <t(t;T(t + 0,s + 0)<p,if)dT. 

By the uniqueness of the solution of the weak form of our evolution equation, we have 
T(t + 6,s + 0)<p = T(t, s)<p, for <p 6 H. 
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Under the periodicity Hypothesis 6.1, the continuity assumptions and the uni- 
form boundedness assumptions of the control system need only to be verified in the 
bounded interval [0,0]. The above lemma shows that the periodicity of the evolution 
operator is given by the periodicity of the corresponding sesquilinear form a. The fol- 
lowing theorem plays a very important role in the study of periodic systems. We give 
its proof in order to remind the reader of the dependency of certain bounds involved. 

Theorem 6.1 ([HI], [H2]) Let be a 9-periodic evolution operator. Then T(-,-) 

is uniformly exponentially stable if and only if there exist an integer n and a constant 
(3 < 1 such that 

(6.1) l|r(n«,0)|| t ( H ) < /?■ 

Proof: a) Let T(-,-) be uniformly exponentially stable; that is, there exist constants 
Af,u> > 0 such that 

||T(M)||l ( //) < (fs) G Aoo(O). 

Therefore, if we take n large enough such that Mexp{— u>n&} < 1, and let (3 = 
M exp{— um0}, we have that (6.1) holds. 

b) Suppose (6.1) holds. Let C be a constant such that for all 0 < s < t < n0, 
lir(M)ll < C. Now for any 0 < s < t < oo and t — s > 9, we can find integers k and 
m such that 

k < - — -p < k + 1, (m — 1)0 < s < mO. 

n9 

Therefore, by the semi-group property of T(-,-), we get 

T(t,s) — T(t,(nk + m)9)T((nk + m)0, m0)T(m0, s). 

By definition of k and m, we have m0 — s < nO and t — ( nk + m)0 < n0; then by 

(6.1) , we have 

mM)|| t (H, < CV < C 2 e k '-^ 

Since (} < 1, we have log/? < 0; therefore we find 

||T(M)lk(//) < C 2 e _log/? • e (fc+1)log/3 < 
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with M = C 2 f 0 and u = — log 0/nO. 

Since the evolution operator T is also n0-periodic, therefore we can assume 
without loss of generality that \\T(0, 0)||l(H) < 0 < 1. The following lemma is an 
interesting consequence of the above theorem. 

Lemma 6.2 For a periodic system, if the stabilizability and the detectability assump- 
tions are satisfied, then we can find 6-periodic operator valued functions K(-) and 
\P(-) such that the evolution operator and Tip corresponding to the perturbation of 
the evolution operator T by BK, 'PVL 1 / 2 are also uniformly exponentially stable. 

Proof: Suppose K, <P are operator valued functions in the stabilizability and de- 
tectability assumptions. Let Tk,T p be the evolution operators corresponding to per- 
turbation of T by BK and tPlT 1 / 2 respectively. Without lost of generality, we can 
1 assume that there exists a constant 0 < 1 such that 

(6.2) ||r*(M)|| < 0 \ 

(6.3) ||T*(M)|| < 0. 

Now define ^-periodic operator valued functions K , V P as K(t) = K ( (), »(J) = *(«), 
j for t G [0,0), and extend periodically for t > 6. Then we have 7V(0, 0) = 7/c(0, 0), 

Tvp(0, 0) = 7^(0, 0); therefore (6.2), (6.3) still hold for the new evolution systems. By 
( Theorem 6.1, we conclude that TV, Ty are also uniformly exponentially stable. 

\ As a consequence of this lemma, we can assume without loss of generality that 

} the operator valued functions K , ^ in Hypotheses 4.5 and 4.6 are 0-periodic. In fact 

we can make the following equivalent assumptions: 

> 

Hypothesis 6.2 There exist a constant 0 < 1 and 6-periodic operator valued func- 
tions K, \P, such that if T K , Tip are the evolution operators corresponding to the per- 
turbation of the evolution operator T by BK , iPVL 1 / 2 respectively, then 

\\Tk{@, 0)|| < 0, ||T»(M)|| <0- 
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For the remainder of this section, we shall assume Hypotheses 6.1 and 6.2 
hold and we focus on the uniform stabilizability and the uniform detectability of the 
approximate control systems. Let K N (t) = Pff K(t), ty N (t) — Pjfty and Tff,Ty 
be the evolution operators corresponding to the perturbations of T N by B N K N and 
'P^(VF^) 1 / 2 respectively. As we have seen in the proof of the Theorem 6.1, the 
constant M and decay rate u depend only on the uniform bound of T and constant (3. 
We already know (use the arguments of Theorem 4.4 and uniform boundedness of the 
perturbations) that \\Tff(t, s)||l(t/ n ) an d ||7# r (M)||L(n*) can be uniformly bounded 
by a constant C independent of N for all 0 < s < t < 0; therefore to obtain uniform 
stabilizability and uniform detectability, we only have to show (see b) of the proof of 
Theorem 6.1) that we can find ft < 1 and N 0 such that for all N > iV 0 , we have 

(6.4) lir A A '(«,o)|| t( „», < k 

(6.5) l|r£'(0,o)|| I . ( „» l < A 

If, on the other hand, the convergence of T^(t,s)Pfl to Tx(t,s) (T§ (t,s)Pfj to 
Ty(t,s)) is in the operator norm, then we can readily establish that (6.4), (6.5) hold 
for N large enough. The following theorem is very useful in the proof of this desired 
convergence. 

Theorem 6.2 Let //, V be Hilbert spaces as defined in section 4 with V compactly 
embedded in H. Consider a sequence of bounded linear operators T N defined on H 
and bounded linear operator T defined on H. Suppose the range ofT N and T are in 
V , and the following conditions hold: 

(i) There exists a constant C such that 

\\T N \\l(h,v) < C, ||T|| t ( H( v) < C\ 

(ii) For any y € H , we have 

II T N <? - 7V||„ 0, UT-'V - TVIIh ^ 0, 

as N — * oo. 
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Then the convergence of the sequence of operators T N to T is in the operator norm, 
that is || T N — T\\l(h) — ► 0, as N — ► oo. 

Before we give the proof of this theorem, let us state a useful lemma. 

Lemma 6.3 Consider a nonincreasing sequence of compact sets Ek, Ek 2 Ek+i, 

k = 1,2,3, If we have fl^i Ek = {0}, then for each e > 0, we can find k Q large 

enough such that for all k > ko, Ek is a subset of a ball B( 0, e) in H defined by 
B{0,e) = {' P eH\M„<e}. 

Proof: Suppose there exists e > 0 such that for every k, we can find <p>k € Ek, such 
that ||<£>fc||f/ > e. Since the sequence {<*?&} is in E\ which is compact, we can assume 
that tpk converges to an element <p in E\. Obviously ||^>|| h > e/2. On the other hand, 
<p must be in Ek for all k , therefore <p is also in the set E = HfeLi Ek- But since 
E = {0}, this is a contradiction. 

Proof of Theorem 6.2: i) By definition of the operator norm, we have 
II T N - T|| 1( „) = sup {\\T N f - 7>|Ih I f € H, I^IIh < l) . 

Now let us define the set Fk as 

OO 

U {T N V - 7V | € //, Mj/ < 1}. 

N=k 

Let Ek be the H closure of Fk- The sequence of operators T N converges to T in 
operator norm if and only if for all t > 0, we can find k 0 such that for all k > ko, we 
have Ek Q B( 0, e). 

ii) We observe that Ek is a closed set in H and hence in V, and by our 
assumption Ek C V is bounded in V norm, in fact 

Ek cW\<pev,Mv<2C}. 

Therefore, Ek is a compact set. By definition of Ek, we know that {Ek} is a nonin- 
creasing sequence of compact sets. 
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iii) Let us define E = p|fcLi -E*. Suppose <p 6 E] since yp is in the closure of Fk 
for all k , we can find a sequence -) N with ||7 ^||h < 1, such that yp N = (' T N — T)^ N 
converges to yp. Therefore, for any € H , we have 

<<?,*/> >H = lim < <p N , >H 

N—> oo 

= lim < 7 N ,(T* N -T*)tl>>„. 
yv— t-oo 

Since 7 ^ is uniformly bounded in N and ( T* N — T*)ip goes to zero as N — » 00 , we 
have < yp, >h= 0 for all € H, and therefore yp = 0 . 

Using the previous lemma and i), we obtain HT^ — T\\l(h) — * ► 0, as N — + 00 . 

Now recall the definition of the approximate control systems defined in section 
4, and consider T N = Tj^(0,0)P ^ , T = Tk(6, 0). By the results of the section 4, 
we can easily verify that the assumptions of Theorem 6.2 are satisfied. Therefore, we 
have 

\\TZ(e,0)Pg-T(0Mnm^<>, 

as N — * 00 . Now let /3 be the constant in Hypothesis 6.2; letting t = 1 — /?, we can 
find No large enough such that for all N > N 0 we have 

||r«(o,o)p^-r(«,o)|U (H) <^. 

Therefore, for all N > No, we have 

11^(0, o)IU(*») < l|r K («,o)||i( H) + ||rj5'(9,o)p# - iw(«,o)|| 1{H) 

e 

“ 1 ~~ 2 ‘ 

We summarize our discussion in the following theorem: 

Theorem 6.3 Lei H, V be the Hilbert spaces used in the section 4 and assume that V 
is compactly embedded in H . Suppose that Hypotheses 4-l~4-3, ^.7 and 6.1, 6.2 hold. 
Let H N C V be the finite dimensional approximation spaces and let the approximate' 
control system be defined as in section 4 ■ Then we can find No large enough such 
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: that for all N > No, the approximate control systems are uniformly stabilizable and 

i uniformly detectable. v4s a consequence, if Q N , Q are the unique solutions of the 

| Riccati integral equations on the infinite time interval in H N , H respectively, and 

I S N ,S are the evolution operators corresponding to the perturbations of T N , T by 
—B^(R n )~ 1 B* ! ^Q n , —BR~ l B*Q respectively, then we have: 

\\Q rf {t)Pg'p - Q{t)<p\\„ -> 0, II S s (t,s)Pg - S(t,s) \\ L{H) -> 0, 

as N — > oo. The convergences are uniform in [0,0] and for ( t,s ) 6 A(O,0) respec- 
tively. 

We remark that an autonomous system is a particular case of a periodic system; 
therefore the approach used here can also be applied to time invariant systems as 
1 considered in [BK 1] . In the case of parabolic systems with strict V-ellipticity where 

V is compactly embedded in H , the arguments in this section offer an alternative (and 
more succinct) approach to uniform stabilizability/detectability from the dissipativity 
approach of section 5. 

\ 

7 Parabolic partial differential equation control systems: 
An example 

In this section we consider control systems governed by second order parabolic partial 
differential equations with distributed scalar control. We indicate briefly how one 
i formulates the associated control and approximate problems in the framework of 

section 4. For this class of systems we show that one can, under standard assumptions, 
f readily verify the conditions for continuity, ellipticity, stabilizability and detectability 

l required in section 4. For Galerkin type approximation schemes based on spline 

subspaces, Conditions US and UD are readily established. 

Let fl be a bounded open subset of with an infinitely differentiable bound- 
ary r given by a variety of dimension n — 1 and consider the following homogeneous 
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second order parabolic partial differential equation ([LI], [L2, pp. 100]): 

(7- 1 ) 

+c(U)*(U), t>o,(en, 

where £ = (£1, • • • ,£ n ) G 3£ n . Generic boundary conditions are given by 

(7.2) 7(0 E + = *>o,£er, 

i,j= i 

where f/(£) = (t;^), ’ ' ’ , Vn{0) ls the outward unit normal vector at a point £ on the 
boundary T of 0. Note that for all £ 6 I\ if 7(£) ^ 0, we can divide (7.2) by 7(£); 
therefore, we can assume without loss of generality that 7 takes only values 0 or 1. 
We choose as our state space H = L 2 (D); the appropriate choice for V depends on 
the boundary conditions and we consider several special cases. 

(i) Consider the case 7(£) = 0, /?(t,£) = 1, and thus equation (7.2) specifies the 
usual Dirichlet boundary condition. We then define V = Hq(Q,), and a sesquilinear 
form <j\ by 

ij=i a st 

- 1 £ {ut.o^viom + cfuMomj <%■ 

(ii) If we consider the case for 7(£) = 1, /3(t,£) = 0, we obtain a Neumann 
boundary condition. We then choose V = (f l) and note that the integrals in the 
definition of a\ above are also defined for any functions <^>, rp in H 1 ( Q). Therefore the 
sesquilinear form cr 2 for Neumann boundary conditions can be taken as the same as 
for Dirichlet conditions, cr 2 = <T\, and thus only the spaces V are changed. 

(iii) Consider the case 7(£) = 1 and 0(t,() ^ 0 which results in Robin or 
mixed boundary conditions. We again choose V = H 1 ( fl), and define a sesquilinear 
form a 3 by 

= <r 2 {t]V,'l’) + J r P(t,£ MO^iOdC 
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By taking 7(£) = 0 on a part Tj of the boundary and y(£) = 1 on T — Tj, 
we can obtain other mixed boundary conditions. The choice of space V should also 
be modified accordingly. In all the cases above, let V be the appropriate choice of 
Sobolev space (either Hq or H 1 ), and let a denote the corresponding sesquilinear 
form defined on V x V. Then a solution z(t) of equations (7.1) and (7.2) satisfies 

(7.3) < z(t), xp > H = -cr(<; z(t), xp), for all xp G V. 

As usual (7.3) is called the weak form of equations (7.1) and (7.2); see equations (4.7), 
(4.8) and (4.10), (4.11) of section 4. 

The continuity and the ellipticity conditions for the sesquilinear forms <r ,• can 
be characterized by properties of the coefficients a,j, 6,-, c and /?. The standard as- 
sumptions ([L2], pp. 100) for F-continuity and Holder continuity of the sesquilin- 
ear forms cr,- are as follows: For each fixed t, the functions a,j(t, •), &,•(<,•), c(t,-) 
are elements of L°°( fl), while the function /3(t, •) is an element of L°°( T). Further- 
more, for each bounded interval [a, 6], there exist constants C > 0 and 0 < 7 < 1, 
such that each of the coefficients a, t j,&,-,c, /? satisfy the bounds ||/(t, •)||£,°° < C and 
ll/(V) — /(-SjOIIl 00 < C\ t — s|° for t,s in [a, 6], where L°° is L°°(Q ) or £°°(r) as 
appropriate. 

It is easily seen that under these assumptions, the F-continuity and F-Holder 
continuity Hypotheses 4.1 and 4.3 hold for G\ and <r 2 defined above. In the case of 
<73, a boundary integral is involved; but under our assumptions on the smoothness of 
the boundary, the following estimates hold (see ([L2], pp. 17, Theorem 3.2, pp. 23): 
For any ip,xp £ H 1 ( fl), the restrictions of tp,xp to the boundary T belong to L 2 (T). 
Furthermore, there exists a constant C such that 

J^tptpd^ < |MlL 2 (r)IM|L 2 (r) < 

Furthermore, for all e > 0, there exists constant C(e) such that 

IMIl 2 (r) ^ + C^OIMIl 2 ^)- 


i 
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With these estimates, it is readily argued that cr 3 also satisfies the Hypotheses 4.1 
and 4.3 of section 4. 

To assure V -ellipticity we again make standard assumptions: for any bounded 
interval [a, 6], there exists a constant v > 0 such that for all t E [a, b] and £ 6 fl, 

^ £ c, 2 < J2 

«=i «.i= i 

for all ( = (Ci, • • • , Cn) £ Under this assumption, it is readily seen that for each 
bounded time interval [a, 6], there exist constants c 2 > 0 and m, such that 

Re <r,(f;y?,y>) > c 2 ||^||y - m||v?||^, 

for all (f 6 V and t 6 [a, 6]. That is, each of the sesquilinear forms cr,, i = 1,2,3, 
defined above satisfies Hypothesis 4.2. 

In the remaining part of the this section, let H, V be the spaces of functions 
appropriate for a specific problem, and let a be the sesquilinear form defined for that 
problem as above. For the control space U, we choose U = T 2 (fl), with the control 
system being defined by (see (4.10)-(4.12)) 

(7.4) ^- t < >n= -a{t]z(t),rp)+ < B(t)u(t),i b > H , 

for all V’ in V. We choose the cost functional given by (4.12). Here we define the 
operators B(t), W(t), R(t) by the following 

[B{t)v ] (0 = b(t,£)v((), for v E L 2 (fl), 

[W(t)v] (£) = u>(t, £)v>(0> for V 5 e T 2 (H), 

[R(t)u] (0 = r(t, 0»(0> for v e ^ 2 (^)> 

where 6, w, r are scalar valued functions on [t 0 , oo) x fl. In this case, the uniform 
boundedness and the positivity of the corresponding operators can be readily char- 
acterized by conditions on the functions 6, w and r. We assume that for each fixed 
t, b(t, •) and r{i ,•) are elements of Z,°°(fl). As L°°(Q ) valued functions of t , 


52 



these functions are assumed continuous. Furthermore, assume there exists a constant 
C such that for all t E [to? 00 )? the functions b(t,-),w(t, -),r(t,-) satisfy the bound 
||/(t, -)||£,oo(Q) < C. The functions in, r are assumed to be nonnegative and, indeed 
we assume there exists a constant r 0 > 0 such that r(t,£) > r 0 , a.e. in 0, for t > t 0 . 

Under these assumptions, the operator valued functions B,W,R satisfy the 
standing assumptions of sections 2 and 4. It remains to consider Hypotheses 4.5 and 
4.6 (stabilizability and detectability of the original system) as well as Conditions US 
and UD once we have introduced approximations. For the problems considered here, 
we can use the definition of the sesquilinear forms to give sufficient conditions for 
Hypotheses 5.1 and 5.2 (and hence Hypotheses 4.5 and 4.6) to hold. To this end, 
we assume that there exist constants p > 0 and p > 0 such that for t E [to, oo), 
> /b KUOI > /b a.e. in fo- 
under this assumption, it is readily seen that there exists constants / > 0 and 
k > 0, such that each of the sesquilinear forms a, satisfies for <p E V 

<r(t i V 5 , V 5 ) + 1 < B(t)B*(t)tp, <p>h> fc|MI h> 

and 

° {t;<p,<p) + l <W(t)(p,ip> H > k\\(p\\ 2 H . 

Thus, Hypotheses 5.1 and 5.2 hold with K(t) = lB*(t) and 'F(t) = 

For our approximate systems, we choose approximation spaces H N and U N 
as in sections 3 and 4 generated by finite element or spline basis elements chosen so 
that H n C V and U N C U yield the desired convergence properties for P$ and Py 
respectively (see [C, Chaps. 2, 3], [B], [Sc]). The approximating systems are then 
defined as in section 5. It follows immediately that Hypothesis 5.3 holds and hence 
the conclusions of Theorem 5.2 are valid for the class of examples considered in this 
section. 

We note that under periodicity assumptions we could have applied the alter- 
native approach of section 6 to these examples since (see [A]) both V = /f 1 (H) and 
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V = Hq(SI) embed compactly in H = L 2 (il). 

In some of our related efforts, we have numerically tested the ideas presented 
in this paper on one dimensional versions of the example of this section. In these 
examples = (0, 1) and we have to date used either linear or cubic B-splines to gen- 
erate the approximation spaces H N and U N . (In fact, when fl is a parallelepiped, the 
above theory still is applicable and tensor products of one dimensional elements are a 
good choice for approximation elements.) We have considered several examples with 
time dependent periodic coefficients; for these examples we could use eigendirection 
analysis (see [W]) to give an analytic analysis for the feedback control problems. The 
resulting analytical solutions were used for comparison with the numerical solutions 
obtained using software implementations based on the theory developed in this paper. 
Quite satisfactory results were obtained and, as noted in the Introduction, these are 
being detailed in a separate manuscript under preparation. 

In concluding we note that the theory in this paper is also applicable to higher 
order parabolic systems (as well as to some boundary damped hyperbolic systems 
[BKS], [BKW]). In particular one dimensional Euler-Bernoulli beam models with 
Kelvin- Voigt damping satisfy (see [BIl]) the strong ellipticity assumptions needed in 
the theory developed above. While boundary control (as treated in [BI2]) for such 
models constitute an obvious class of problems, distributed control as treated in this 
paper is essential in cases where nonuniform piezoelectric layers along the beam are 
used to implement the feedback controls. 
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